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f^*) , The low energy effective action of gravity in any even dimension generally acquires 

O I non-local terms associated with the trace anomaly, generated by the quantum fluctuations 

of massless fields. The local auxiliary field description of this effective action in four 
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dimensions requires two additional scalar fields, not contained in classical general relativity, 
which remain relevant at macroscopic distance scales. The auxiliary scalar fields depend 
upon boundary conditions for their complete specification, and therefore carry global 
information about the geometry and macroscopic quantum state of the gravitational field. 
The scalar potentials also provide coordinate invariant order parameters describing the 
conformal behavior and divergences of the stress tensor on event horizons. We compute the 
stress tensor due to the anomaly in terms of its auxiliary scalar potentials in a number of 
concrete examples, including the Rindler wedge, the Schwarzschild geometry, and de Sitter 
spacetime. In all of these cases, a small number of classical order parameters completely 
determine the divergent behaviors allowed on the horizon, and yield qualitatively correct 
global approximations to the renormalized expectation value of the quantum stress tensor. 
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I. INTRODUCTION 

Gravitation is a macroscopic phenomenon, governing the behavior of matter at scales from a few 
cm. up to the largest scales of cosmology, 10^^ cm. On the other hand quantum phenomena are 
associated generally with microscopic distances, on the order of atomic scales, 10^^ cm. or smaller. 
The physical and mathematical frameworks used to describe physics at these disparate scales are 
quite different as well, and at first sight appear even to be irreconcilable. In general relativity 
gravitation is founded on the Equivalence Principle, which finds its mathematical expression in 
differential equations which transform covariantly under a local change of frame at every spacetime 
point. In classical physics this idealization of spacetime as a local differentiable structure is quite 
natural. Yet because of quantum fluctuations it is very likely that the classical description of 
spacetime as a smooth pseudo-Riemannian manifold fails at very short distance scales, of the order 



of the Planck length, Lpi = ^JhGj(? ~ 10^^^ cm. 

Perhaps less widely appreciated is that quantum theory also raises questions about the purely 
local, classical description of the gravitational interactions of matter, even at macroscopic distances. 
This is because of the peculiarly quantum phenomena of phase coherence and entanglement, which 
may be present on any scale, given the right conditions. Indeed macroscopic quantum states are 
encountered in virtually all branches of physics on a very wide variety of scales. Some of the better 
known examples occur in non-relativistic many-body and condensed matter systems, such as Bose- 
Einstein condensation, superfluidity, and superconductivity. These quantum coherence effects due 
to the wave-like properties of matter have their close analogs in relativistic field theories as well, in 
spontaneous symmetry breaking by the Higgs mechanism in electroweak theory, and chiral quark 
pair condensation in low energy QCD. In each of these cases the ground state or "vacuum" is a 
macroscopic quantum state, described at very long wavelengths by a non-vanishing quasi-classical 
order parameter, in a low energy effective field theory. 

The strictly classical theory of general relativity would seem to preclude any incorporation of 
macroscopic coherence and entanglement effects of quantum matter, on any scale. Taking account 
of the quantum wave-like properties of matter, and its propensity to form phase correlated states 
over macroscopic distance scales at low enough temperatures and/or high enough densities requires 
at the least a semi-classical treatment of the effective stress-energy tensor source for Einstein's 
equations. The corresponding one-loop effective action for gravity can be determined by the meth- 
ods of effective field theory (EFT). Our principal purpose in this paper is to demonstrate that 



non-local macroscopic coherence effects are contained in the low energy EFT of gravity, provided 
that Einstein's classical theory is supplemented by the contributions of the one-loop quantum trace 
anomaly of massless fields. 

The low energy EFT of gravity is determined by the same general principles as in other contexts 
m , namely by an expansion in powers of derivatives of local terms consistent with symmetry. Short 
distance effects are parameterized by the coefficients of local operators in the effective action, with 
higher order terms suppressed by inverse powers of an ultraviolet cutoff scale M. The effective 
theory need not be renormalizable, as indeed Einstein's theory is not, but is expected nonetheless 
to be quite insensitive to the details of the underlying microscopic degrees of freedom, because 
of decoupling [l|. It is the decoupling of short distance degrees of freedom from the macroscopic 
physics that makes EFT techniques so widely applicable, and which we assume applies also to 
gravity. 

As a covariant metric theory with a symmetry dictated by the Equivalence Principle, namely 
invariance under general coordinate transformations, general relativity may be regarded as just 
such a local EFT, truncated at second order in derivatives of the metric field gab{x) ,2]. When 
quantum matter is considered, the stress tensor becomes an operator whose expectation value acts 
as a source for the Einstein equations in the semi-classical limit. Fourth order terms in derivatives 
of the metric are necessary to absorb divergences in the expectation value of the stress tensor in 
curved spacetime, but the effects of such higher derivative local terms in the gravitational effective 
action are suppressed at distance scales L ^ Lpi in the low energy EFT limit. Hence surveying 
only local curvature terms, it is often tacitly assumed that Einstein's theory contains all the low 
energy macroscopic degrees of freedom of gravity, and that general relativity cannot be modified at 
macroscopic distance scales, much greater than Lpi, without violating general coordinate invariance 
and/or EFT principles. However, this presumption should be re-examined in the presence of 
quantum anomalies. 

When a classical symmetry is broken by a quantum anomaly, the naive decoupling of short and 
long distance physics assumed by an expansion in local operators with ascending inverse powers 
of M fails. In this situation even the low energy symmetries of the effective theory are changed 
by the presence of the anomaly, and some remnant of the ultraviolet physics survives in the low 
energy description. An anomaly can have significant effects in the low energy EFT because it is 
not suppressed by any large energy cutoff scale, surviving even in the limit M ^ co. Any explicit 
breaking of the symmetry in the classical Lagrangian serves only to mask the effects of the anomaly, 



but in the right circumstances the effects of the non-local anomaly may still dominate the local 
terms. A well known example is the chiral anomaly in QCD with massless quarks, whose effects are 
unsuppressed by any inverse power of the EFT ultraviolet cutoff scale M ~ Aqcd- Although the 
quark masses are non-zero, and chiral symmetry is only approximate in nature, the chiral anomaly 
gives the ^dominant contribution to the low energy decay amplitude of vr*^ — > 27 in the standard 
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model 



a contribution that is missed entirely by a local EFT expansion in pion fields. Instead 



the existence of the chiral anomaly requires the explicit addition to the local effective action of a 



m. 



non-local term in four physical dimensions to account for its effects 

Although when an anomaly is present, naive decoupling between the short and long distance 
degrees of freedom fails, it does so in a well-defined way, with a coefficient that depends only on the 
quantum numbers of the underlying microscopic theory. In fact, since the chiral anomaly depends 
on the color charge assignments of the short distance quark degrees of freedom, the measured 



ymg microscopic 



low energy decay width of vr^ -^ 27 affords a clean, non-trivial test of the under 
quantum theory of QCD with three colors of fractionally charged quarks [l|, U, la]. The bridge 
between short and long distance physics which anomalies provide is the basis for the anomaly 
matching conditions [2|]. 

In curved space an anomaly closely related to the chiral anomaly also appears in massless 
quantum field theory a 1^1 • This conformal or trace anomaly provides us with additional infrared 
relevant terms that do not decouple in the limit Mp^ = Lpi -^ 0, and which should be added to 
the Einstein-Hilbert action of classical relativity, to complete the EFT of low energy gravity. 

That the trace anomaly terms are necessary for the low energy completion of the EFT of 
gravity may be seen from the classification of possible terms in the gravitational effective action 
according to their behavior under global Weyl rescalings 10]. The terms in the classical Einstein- 
Hilbert action scale with positive powers (~ L"^ and ~ L^) under rescaling of distance, and are 
clearly relevant operators of the low energy description. The non-local anomalous terms scale 
logarithmically (~ logL) with distance under Weyl rescalings, and therefore should be considered 
marginally relevant. Unlike local higher derivative terms in the effective action, which are either 
neutral or scale with negative powers of L under global Weyl rescalings of the metric, the anomalous 
terms cannot be discarded in the low energy, large distance limit. The addition of the anomaly 
term(s) to the low energy effective action of gravity amounts to a non-trivial infrared modification 
of general relativity, fully consistent with both quantum theory and the Equivalence Principle jl^] . 

The anomalous terms in the effective action lead to corresponding additional terms in the stress 



tensor and equations of motion of the low energy EFT of gravity. These terms are most conveniently 
expressed in terms of two auxiliary scalar fields which permit the non-local effective action of the 
trace anomaly and its variations to be cast into local form. The additional scalar degrees of freedom 
in these terms take account of macroscopic effects of quantum matter in gravitational fields, which 
are not contained in the purely classical, local metric description of Einstein's theory. 

In this paper we give explicit formulae for the stress tensor due to the trace anomaly of massless 
fields to be added to the classical Einstein equations in the low energy EFT limit, and expose its 
macroscopic effects in a number of familiar geometries. The most interesting effects are associated 
with geometries that have horizons, boundaries, or non-trivial global topologies. For example, the 
two new auxiliary fields in the effective action generically diverge as the apparent horizon of Rindler, 
Schwarzschild, and de Sitter spacetimes is approached. This behavior of the scalar auxiliary fields 
gives a coordinate invariant semi-classical order parameter description of possible divergences of 
stress tensors of quantum field theories on spacetime horizons, enabling an understanding of such 
divergences in terms of the global properties of the geometry. The Casimir energy density between 
two conducting plates is an example of a bounded space without horizons or curvature which also 
admits a simple description in terms of the scalar auxiliary fields. The various examples taken 
together suggest that the degrees of freedom contained in the scalar potentials induced by the trace 
anomaly should be regarded as semi-classical macroscopic order parameter fields (condensates), 
whose non-vanishing values are connected with non-trivial boundary conditions, horizons, or the 
topology of spacetime. 

The paper is organized as follows. In the next section we review the auxiliary field form of 
the effective action of the trace anomaly in two dimensions, recovering results for the stress tensor 
of conformal fields in the two dimensional Rindler wedge, Schwarzschild and de Sitter metrics. 
In Section III we discuss the general form of the effective action in four dimensions, and give 
the additional terms in the stress tensor in terms of the auxiliary fields, as well as the conserved 
Noether current associated with it. In Section IV we consider the effects of these additional terms 
and the new local scalar degrees of freedom of low energy gravity they represent in flat, Rindler 
and the conformally flat de Sitter spacetime. In Section V we discuss the approximate conformal 
symmetry which pertains in the vicinity of any static, spherically symmetric Killing horizon and 
give the explicit form of the auxiliary fields in the Schwarzschild case. We characterize the possible 
singularities of the stress tensor on the Schwarzschild horizon, and find approximate stress tensors 
for the Hartle-Hawking, Boulware, and Unruh states, comparing our results with numerical results 
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for {T^) of massless fields of various spins. We conclude in Section VI with a summary of our 
results and their implications for quantum effects in gravity. The Appendix contains the explicit 
formulae for the components of the stress-energy in static spherically symmetric spacetimes, upon 
which the Schwarzschild and de Sitter results are based. 

II. EFFECTIVE ACTION AND STRESS TENSOR IN TWO DIMENSIONS 

The effective action for the trace anomaly in any even dimension is non-local when expressed in 
terms only of the metric Qab- It can be made local by the introduction of scalar auxiliary field(s). 
To illustrate this procedure and introduce the general framework which we will use extensively in 
four dimensions, we consider first the somewhat simpler case of two physical dimensions. In d = 2 
the trace anomaly takes the simple form [9|, 

(7^a") = ^^, {d = 2) (2.1) 

where A^ = Ng + Np is the total number of massless fields, either scalar (A5) or fermionic (Np). 
The fact that the anomalous trace is independent of the quantum state of the matter field(s), and 
dependent only on the geometry through the local Ricci scalar R suggests that it should be regarded 
as a geometric effect. However, no local coordinate invariant action exists whose metric variation 
leads to (|TT|) . 

A non-local action corresponding to (|2.1|) can be found by introducing the conformal parame- 
terization of the metric, 

gab = e'^^gah , (2.2) 

and noticing that the scalar curvature densities of the two metrics gab and gab are related by 

Ryf^ = Ryf^-2yf^Ua, {d = 2) (2.3) 

a linear relation in a in two (and only two) dimensions. Multiplying 1)2. 1() by ^/—g, using (|2.3|) 
and noting that y/—g{T^'^) defines the conformal variation, ST^"^' /5a of an effective action r*^^\ we 

conclude that the a dependence of T^"^' can be at most quadratic in a. Hence the Wess-Zumino 

(2) 
effective action in two dimensions, T^^ is 

T^Sz[9-^<^\ = ^jd^^V^{-aUa + Ra). (2.4) 



This action functional of the base metric gab and the Weyl shift parameter a may be regarded as 
a one- form representative of the cohomology of the local Weyl group in two dimensions [li| . This 
means that r^y^[^;(7] is closed under the co-boundary (i.e. anti-symmetrized) composition of Weyl 
shifts, 

A., o T^^^[g; a,] = T^^ziae^'"' ! ^2] - rij^^^e^'^^; ^1] + rg.^^; ^1] - rjjy^; ^2] = o . (2.5) 

(2) 

but is non-exact, in the sense that ^yyz cannot itself be written as the co-boundary shift of a local 



action functional. The relation (|2.5|) following directly from the hermiticity of D is exactly the 
con: 

.(2) 



.(2) 



irectly 

m 



Wess-Zumino consistency condition for r^^[5r;o"] 

(2) 
Although r^^ cannot be written as a coboundary shift A^- of a local single-valued covariant 

scalar functional of the metric (jab, it is straightforward to find a non-local scalar functional Sanom[g] 

such that 

r'^z [g;a] = A^o 5£L M ^ S^^^ [g] - S^l^ [g] . (2.6) 

Indeed by solving (|2.3|) formally for o", and using the fact that y/—g\I\ = y/—g\I\ is conformally 
invariant in two (and only two) dimensions, we find that TyJ^ can be written as a Weyl shift (|2.6() 
with 

S^LM = ^I d^^V^I d'x' V^R{x) n-\x,x') R{x') , (2.7) 

and \Ii~^{x,x') denoting the Green's function inverse of the scalar differential operator D. The 
parameter Q^ is —N/6 if only matter fields in a fixed spacetime metric are considered. It becomes 
(25 — A^)/6 if account is taken of the contributions of the metric fluctuations themselves in addition 

n 

to those of the N matter fields, thus effectively replacing A^ by A^ — 25 [12|- In the general case, 
the coefficient Q^ is arbitrary, and can be treated as simply an additional free parameter of the low 
energy effective action, to be fixed by experiment. 

The anomalous effective action (|2.7j) is a scalar under coordinate transformations and there- 
fore fully covariant and geometric in character. However since it involves the Green's function 
n^^(x,x'), which requires boundary conditions for its unique specification, it is quite non-local, 
and dependent upon more than just the local curvature invariants of spacetime. The non-local and 

(2) 

non-single valued functional of the metric, Sanom may be expressed in a local form by the standard 
method of introducing auxiliary fields. In the case of 1)2. 7|) a single scalar auxiliary field, ip satisfying 

-Dip = R (2.8) 



is sufficient. Indeed, varying 

S'^^^omW.V] ^ ^y d^xV^ (Va(^V>-2i?(^) (2.9) 

witli respect to ip gives tfie Eq. of motion (|2.8|) for tlie auxiliary field, which when solved (formally) 
hy if = —\I\~^R and substituted back into Sanom.[g] 'f] returns the non-local form of the anomalous 
action H2.7() . up to a surface term. The non-local information in addition to the local geometry 
which was previously contained in the specification of the Green's function n^^(x,x') now resides 
in the auxiliary local field (/^(x), and the freedom to add to it homogeneous solutions of H2.8() . 
The variation of ()2.9|) with respect to the metric yields a stress-energy tensor. 



(2'] 

1 Sg" 



^afe b;v'] = -^= — T-^ 



^ V7 V7 ^. , „ ^,. -.VT ^.^.V7 ,.^ , -„ .V7 ^.^.v7c,.^ (2.10) 



-VaVfe99 + gabn^p- -(Va(/?)(Vfe(/?) + -gab (Vc(/?)(VV) 
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which is covariantly conserved, by use of (|2.8jl and the vanishing of the Einstein tensor, G, 
Rab — Rgab/"^ = in two (and only two) dimensions. The classical trace of the stress tensor, 



g'^^TS^[g;^]=^U^ = -^R (2.11) 



reproduces the quantum trace anomaly in a general classical background (with Q^ proportional to 
h). Hence (|2.9() is exactly the local auxiliary field form of the effective action which should be added 
to the action for two dimensional gravity to take the trace anomaly of massless quantum fields into 
account. 

The full effective action of two dimensional gravity is 



(2) 

where S^i is the local classical Einstein-Hilbert action, 

(2)r„i _., / ^2^,^— 7? , \ / ^2 



5g) b,(/.] =5£U[9,y.]+4')b], (2.12) 



S'ci' [g]=l j d'x,/^R + \j d'x./^, (2.13) 

which are the only terms we would have written down according to the usual EFT approach, 
expanding in strictly local terms up to and including dimension two which transform as scalars 
under general coordinate transformations. Under global rescalings of the metric ()2.2|) with a = ao, 
a spacetime constant, the volume term scales like e^'^" and its coefficient A has positive mass 
dimension two. Hence it is clearly a relevant term in the effective action at large distances. The 
integral of the Ricci scalar is independent of o"o, with a coefficient 7 that is dimensionless, while 



the integrals of all higher local curvature invariants scale with negative powers of e*^", multiplied 



by coefficients having negative mass dimensions. These terms are neglected in (|2.13|) since they are 
strictly irrelevant in the low energy EFT limit. 

The anomalous term (|2.7|) . or (|2.9|) scales linearly with (Tq, i.e. logarithmically with distance, 
and is therefore a marginally relevant operator in the infrared. Even if it is not included at the 
classical level it will be generated by the one-loop effects of massless fields, which do not decouple 
at any scale. Hence it should be retained in the full low energy effective action 1)2. 12() of two 
dimensional gravity. Moreover, since the integral of i? is a topological invariant in two dimensions, 
the classical action 1)2. 13() contains no propagating degrees of freedom whatsoever, and it is Sanom 
which contains the only kinetic terms of the low energy EFT. In the local auxiliary field form 
1)2. 9|) . it is clear that Sanom describes an additional scalar degree of freedom cp, not contained in 

(2) 

the classical action S^^ . This is refiected also in the shift of the central charge from N — 26, which 
would be expected from the contribution of conformal matter plus ghosts by one unit to A^ — 25. 

Extensive study of the stress tensor and its correlators arising from this effective action estab- 
lished that the two dimensional trace anomaly gives rise to a modification or gravitational "dressing" 
of critical exponents in conformal field theories at second order critical points [1^. Since critical 
exponents in a second order phase transition depend only upon fiuctuations at the largest allowed 
infrared scale, this dressing is clearly an infrared effect, independent of any ultraviolet cutoff. These 
dressed exponents are evidence of the infrared fiuctuations of the additional scalar degree of free- 
dom (p which are quite absent in the classical action (|2.13jl . The appearance of the gravitational 
dressing exponents and the anomalous effective action 1)2. 7|) itself have been confirmed in the large 
volume scalin g lim it of two dimensional simplicial lattice simulations in the dynamical triangulation 
approach [ij, |lJ]. 

The formal similarity between (|2.4j) and (|2.9j) suggests that the introduction of the local auxiliary 
field if has simply undone the steps leading from (|2.4() to ()2.7() . which eliminated the conformal 
factor a. By comparing (|2.3|) and (|2.8|) . we observe that the metric 

gab = e-'^gab (2.14) 

is a metric with zero scalar curvature, R = 0, conformally related to the physical metric gab- Hence 
ip{x) does parameterize the local Weyl transformations of the metric along the same fiber as cr(x) 
in the conformal parameterization ()2.2() . The difference between the two representations is that 
whereas gab was regarded as a fixed base metric related to the physical metric gab by 1)2. 2|) . the 
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possibility of adding homogeneous solutions to (|2.8|) shows that (p and hence gab defined by (|2.14|) 
are not unique, and hence no gab plays any privileged role over any other. The only physical metric 
is gab- The non-uniqueness of c^ allows for additional degrees of freedom at boundaries or coordinate 
singularities of gab^ which need not be associated with any smooth gab- The completely covariant 
local form of the action H2.9() makes it clear that these effects, associated with the boundary condi- 
tions on the spacetime scalar ip satisfying ()2.8|) are not coordinate artifacts which can be removed 
by reparameterizations of coordinates in the physical spacetime gab- It is the freedom to specify 
boundary conditions for the classical auxiliary field 99, or equivalently add homogeneous solutions 
to the non-local Green's function n^^(x,x') which implies the possibility of global macroscopic 
effects of the anomaly. 

Both the auxiliary field equation of motion (|2.8|) and stress tensor ()2.1U|) are left invariant by 
the constant shift, 

tp^tp + ifQ. (2.15) 

The underlying reason for this shift symmetry is that the integral, 

X = ^ j d^Xy/^R (2.16) 

is a topological invariant in two dimensions. For Euclidean signature metrics x = Xe is the Euler 
number of the manifold. This implies that in two (and only two) dimensions ^J—gR can be expressed 
as a total derivative, or equivalently, 

R=-^^a{^/^^'') = Va^'', id = 2). (2.17) 

V-9 



The topological current fi", analogous to the Chern-Simons current for the gauge theory chiral 
anomaly J15l.lla|. depends upon the choice of gauge or coordinate frame and is non-unique, although 
its line integral around a closed contour, 

(j> yf^n''nads= (p .f^Q."'eabdx^= y/^RdPx (2.18) 

Jdv Jav Jv 

is a coordinate invariant scalar quantity. Here we have used the notation, na = ^ab-^ for the 

normal to the boundary dV of the two dimensional region V and Cab is the standard alternating 

tensor in two dimensions (with ei2 = +!)• The existence of ri" satisfying 1)2. 17() implies that the 

action (|2.9p may be rewritten in the alternative form, 

Sanom,2[9;v]=^ [ d^X./^{V''ip + 2n'')Vap-^<{ V^^P^^nads, (2.19) 

167r Jy 8tt Jqv 
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leaving (|2.8|) . (|2.1()|) and (|2.11|) unchanged. Since up to the last surface term the anomalous action 
depends on Lp only through its derivatives, there is a frame dependent Noether current corresponding 
to the global Weyl rescaling H2.15() . 

r = V^ip + 17" , (2.20) 

which is covariantly conserved, 

Var = ^daiV^r) = n^ + R = o, (2.21) 

by virtue of dHJ and (|^T7|l . 

Within the framework of the low energy effective action 1)2. 12() all fields may be treated classically, 
with h contained implicitly only in the c-number coefficients 7, A and Q^. In this (semi-)classical 
framework, the auxiliary field c/? is a classical scalar potential which contains information about 
the macroscopic quantum state and topological boundary effects. The geometric information ip{x) 
contains cannot be reduced to local geometric invariants constructed from the curvature tensor, 
its derivatives and contractions at x. Instead the massless correlation function \3^^{x,x'), which 
grows logarithmically with the invariant distance between the spacetime points x and x', and the 
existence of a conserved global charge corresponding to the Noether current J" imply that ip is a 
long range field depending upon the global properties of spacetime. It is this long range behavior 

of ip that enables the effective action (|2.12|) to incorporate macroscopic quantum coherence effects, 

(2) 
quite absent from the purely local classical geometric action S^^ . Since the effective action contains 

(2) 
Sanom, these non-local macroscopic quantum coherence effects are incorporated in the low energy 

EFT of two dimensional gravity in a natural way. 

The long range macroscopic effects and the physical meaning of the Noether charge are best 

illustrated by means of a few simple examples. Consider first fiat spacetime in Rindler coordinates, 

ds^ = -dt^ + dx^ = -p^djf + dp^, (2.22) 

with X = pcoshr/, t = psmhr]. The Killing vector dri generates Lorentz boosts and this vector field 
becomes singular on the light cone, p^ = x"^ — t^ = 0. If the state of the system is boost invariant 
then the semi-classical auxiliary field can be assumed to be independent of 77. With ip = 'p{p), (|2.8() 
has the solution (for finite non-zero p), 

(^ = 2gln(-^j , (2.23) 
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with q and pQ arbitrary constants. 

If g 7^ the singularity of y? at p = gives a delta function contribution to the scalar curvature 
at the origin of the Euclidean metric obtained by replacing in H2.22() . r/ -^ iO, i.e., 

nc^ = 47rg^, (2.24) 

P 

so that 1)2. 23() should properly be considered a solution of ()2.8|) only in the distributional sense. If 

the point at p^^ = is included in the manifold then there is also a delta function contribution to 

(|2.24|) at infinity. In the Euclidean signature metric, the operator D is the usual two dimensional 

Laplacian and ()2.24|) is Laplace's equation for the electrostatic potential with a point charge of 

magnitude —q at the origin, and +q at infinity. A non-zero value of q corresponds to a topological 

defect on the Minkowski light cone, with the metric e"'^ ds^ , conformal to the original metric (|2.22|1 

possessing a conical singularity at p = for generic g 7^ 0. 

Substituting (|2.23|1 into the stress tensor H2.1U|1 . we find in the (ry, p) coordinates, 

^'^s^ ("::)■ 

The case q = {) corresponds to the Minkowski vacuum with no singularity on the light cone, t = itx. 
However when g = 1, 

^2- 



e 



-%=, ds' = pI {-drj^ + y-)= Poi-dri' + df) , (2.26) 



with ^ = ln(/9//9o)) is again flat. The vanishing Euler number of this flat metric may be regarded 
as resulting from the cancellation of the — 1 Euler charge at the origin from (|2.8|) , (|2.1()j) and (|2.24|) 
with the +1 Euler number of a circular flat disc with boundary dV at large but finite radius, as 
the boundary is taken to infinity in the original metric. In the conformal metric ()2.26|) in (r/, ^) 
coordinates, for q = 1 both singularities at p = and p = 00 are removed to infinity and we again 



obtain a flat spacetime with no singularities and no boundaries. The original flat spacetime (|2.22|) 
has been conformally mapped to another flat spacetime with a different global topology. Indeed 
the Euclidean signature metric dO"^ + d^'^ with 9 27r-periodic is the metric of i? x S*^, with punctures 
at ^ = ibcxD. The Euler charge of ±1 unit at ^ = ±00 is also a consequence of general theorems 
relating the Euler number to the number of fixed points of the Killing field K = ■^, where the 
metric ()2.22() (or its inverse) becomes singular in the original (p, ??) coordinates jl71 |. 

When q = 2 the metric, e"''' ds"^ becomes equivalent to the original flat metric 1)2. 22() by the 
change of variables p -^ Po/p which turns the spacetime inside out, exchanging the coordinate 
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singularities at p = and p = oo. Since the new metric is equivalent to the original metric up 
to a (singular) coordinate transformation, the stress tensor (|2.25j) again vanishes for q = 2. The 
invariance of the physical stress tensor H2.25() under the transformation q ^ 2 — q is evidence for a 
topological two-fold degenerate vacuum in two dimensional gravity. 

For a single free scalar field Q^ = —1/6 and at the value of g = 1, (|2.25l) gives 

1 f -I 



rp b 



(2.27) 



FR 24:TTp^ \ Q -^ 

in the (r/, p) coordinates. This is exactly the expectation value of the stress tensor of a quantum 
scalar field in the boost invariant Fulling-Rindler state [l3|. The preceding discussion shows that 
the divergence of T„* in the Fulling-Rindler state is associated with a unit topological defect on 
the Minkowski light cone, contained in the auxiliary potential cp, corresponding to a puncture at 
p = (equivalently, p = oo) in the analytically continued Euclidean signature metric. A singular 
conformal transformation 1)2. 23() with (7 = 1 is required to transform the Minkowski vacuum to the 
Fulling-Rindler state. Whereas the Hamiltonian defined with respect to the global Minkowski time 
coordinate t is bounded from below and possesses a well-defined vacuum state, the boost generator 
corresponding to the Killing field dj^ changes sign and is unbounded from below. These global 
properties distinguishing the Fulling-Rindler state from the Minkowski vacuum are refiected in the 
topological defect (|2.24|) and divergence of the stress-energy tensor (|2.27)) on the light cone. 

A second example of macroscopic effects of topology of the trace anomaly is provided by two- 
dimensional static metrics of the form, 

ds^ = -fdt^ + ^ = f{-dt'^ + dr*^), (2.28) 

where / = /(r) and dr* = dr/f. This class of metrics include the two-dimensional Schwarzschild 
metric, with /(r) = 1 — 2M/r and de Sitter metric in static coordinates, with /(r) = 1 — H^r^. The 
Killing field K = dt is timelike for / > 0, and the scalar curvature corresponding to 1)2. 28(1 is R = 



— /". Hence in this case a particular inhomogeneous solution to 1)2. 8(1 is <p = ln{—K°'Ka)'^ = 
The conformally transformed line element e~'^ds^ = —dt^ -|-dr*^ is then the fiat optical metric 
The components of the stress tensor due to the anomaly are given by 

N f (p^r^r* f _ ,1 /.2 , -2^ 



iii- 
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^t' = ^[-j + Yf^,r* - \ {^l* + Vi') + i?} , and (2.29b) 

T^ = ^[p,r^-^-!^^ + \^^A. (2.29c) 
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The off-diagonal ffux component is generally time dependent unless we restrict ourselves to at most 
linear functions of time. With this restriction, the general solution of (|2.8|) leading to a stationary 
stress tensor in the Schwarzschild case is: 

where cq, q and p are constants. As in the Rindler case this is a solution to ()2.8|) in the distributional 
sense, since D ^ contains delta function singularities. Indeed with p = the strength of the delta 
function at the horizon, mapped to a point on the Euclidean section t ^>- it for real r, is 47r(g + 1). 
For q = —1 this horizon singularity is cancelled, although that value leads to a divergence in 1)2. 3U() 
as r — > oo. 

The stress-energy tensor for A^ scalar (or fermion) fields in the Schwarzschild case with / = 
1 - 2M/r is 

N ( I fp'^ + q^ 4M2\ 4M1 






TX-^i^^-'-^). (2.310 



1927rM2 / 
N /pM 
967r/ V 4M2 

This general form of the stress tensor for stationary states in the two-dimensional Schwarzschild 
metric was obtained in ref. Q^, from considerations of conservation, stationarity, and the trace 
anomaly. The identification of the arbitrary constants K and Q in Eqs. (3.4) of that work with 
—pq/192'K and [{p — q^ — l]/3847r respectively of the present article shows that the stress tensor 
obtained from the anomaly, (|2.1()|) with (|2.3()|) coincides with that obtained in [2i|. 

The stress tensor (|2.31|) is generally singular at r = 2M, diverging as f~^ as / — > 0. The 
condition for finiteness on the future horizon is 

|i(r,'-T/+2r,-,l<oo a, .-2M. (2.32) 

Substituting ()2.31() into this condition and expanding about r = 2M gives the condition, 

[q — p) = 1 , (finiteness on future horizon) (2.33) 

Finiteness on the past horizon would imply {q + p)"^ = 1. Thus, finiteness on both horizons implies 
1 p^ + q^ = 1, pq = Q and zero flux. Since (|2.33|) cannot be satisfied with both q and p equal to 
zero, the solution ()2.3U|1 must diverge as either r* or t goes to infinity. Note that this is completely 



15 



unlike the previous Fulling- Rindler case in which a finite (in fact, vanishing) stress tensor on the 
horizon is obtained with 99 = 0, which has both zero charge on the light cone, and no divergence at 
infinity. 

A time independent state with zero flux satisfying (|2.3c{|) is the two-dimensional Hartle-Hawking 
state with either p or q equal to zero, and the other eqial to ±1. The stress tensor for this state is 
given by 



Tt' 



rp r 



T 



= , (2.34b) 

HH 

N ( 1 4M2 , , , ^ 

(2.34c) 



HH 967r/ V4M2 r^ 



The degeneracy of the stress tensor for two different values of p or g is similar to the behavior 
obtained in the Rindler case, (|2.25l) . As the horizon r -^ 2M is approached. 



rp b 



HH 6 




^ ' ' - ■ 2M , (2.35) 



where Th = (SvrM)^^ is the Hawking temperature. Thus the stress tensor approaches that of a 
two dimensional perfect gas at T = Th near the horizon in the Hartle-Hawking state. The same 
value is also obtained in the limit r — > 00, showing that this state has infinite total energy. 

Another state obeying the finiteness condition (|2.32j) on the future horizon is the Unruh state 
p = —q = 1/2. Its time inverse, given hy p = q = —1/2 is finite on the past horizon. These 
states have non-zero fiux T^^ either outwardly or inwardly directed from the horizon extending 
to infinity. The non-zero q in these states reflects singularities, i.e. sources or sinks of the current 



(|2.20;) at both r = and r = 00. 

Inspection of 1)2. 31(1 shows that the only state which has vanishing stress tensor as r — > 00 is 
the Boulware state with p = q = 0, which does not satisfy the finiteness condition (|2.32|) at the 
horizon. Unlike the previous Fulling-Rindler example in which finiteness at both the horizon and 
infinty is achieved by one choice of 99 = 0, in the Schwarzschild metric 1)2. 28() because of its different 
topology, one is forced to choose between finiteness on the horizon and falloff at infinity. In the 
Boulware state with p = q = 0, 

tJ' -^-'-^ I I , r^2M, (2.36) 
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as r ^ 2M and / ^ 0, corresponding to the negative of the stress tensor of a thermal distribution 
of massless particles at the blue-shifted local temperature Tioc = Tn/^fJ ■ 

Comparing ()2.36() with (|2.27|) we observe that the behavior of the stress tensor in the Boulware 
state in two dimensional Schwarzschild spacetime is locally the same as that of the Fulling-Rindler 
state in flat spacetime near the corresponding horizons, and indeed with the correspondence, p <-^ 
4AI^/f, the behavior near the horizons in the two cases may be mapped onto each other. However 
the two situations are quite different globally. The Euler number of flat space is zero, and no 
singular behavior of (p is required on either the light cone t = zizx or infinity. The Fulling-Rindler 
state can be obtained only by requiring a non-zero topological charge g' = 1 in the auxiliary field, 
which leads to if of 1)2. 23(1 diverging at both large and small p. 

On the other hand the scalar curvature R = —f" is non-zero in the Schwarzschild case and the 
Euler number of the corresponding non-singular Euclidean signature metric is +1. This non-trivial 
topology of the background spacetime requires an inhomogeneous solution to (|2.8() , such as ip = In f 
in (|2.80|) . The divergence on the horizon of this particular solution to the inhomogeneous equation 
can be cancelled by homogeneous solutions with non-zero p and/or q = —1 satisfying 1)2. 33() . but only 
at the price of a non-zero energy density, and non-trivial topological charge at infinity. Conversely 
the solution of (|2.3()|) with cq = p = q = which vanishes at infinity and produces a transformed 
metric e~^ds'^ which is flat there, necessarily has a diverging stress tensor on the horizon. Thus 
while particular solutions with particular boundary conditions for the auxiliary field can be chosen, 
the impossibility of producing a solution of (|2.3()j) which is well-behaved at both the horizon and 
infinity is a global property of the Schwarzschild spacetime (|2.28|) . which is quite distinct from flat 
spacetime in Rindler coordinates ()2.22|) . The stress tensor of the anomaly and the auxiliary scalar 
(/? captures these global effects of the macroscopic quantum state in terms of the classical solutions 
of ((2S1). 

Similar conclusions obtain in the two dimensional de Sitter spacetime ()2.28j) with / = 1 — H^r^. 
The general solution of (|2.8j) leading to a stationary stress tensor is 

9? = Co + 2qHr* + 2pHt + In / . (2.37) 
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The components of the stress tensor in the general stationary state labelled by {p, q) are 

Tt' = ^{-j{p' + l'-H'r')+^}, (2.38a) 

~ -' ^^' y + q'- H\^) . (2.38b) 



247r/ 
NH"^ pq 
12^7 



T/* = tl^q, (2.38c) 



The generic stress tensor again diverges on the horizon where / = 1—H'^r^ = 0, despite the finiteness 
of the scalar curvature R = —f" = 2H^ there. The condition that the stress tensor be finite on the 
future horizon is again {q — p)'^ = 1. The divergence of the stress tensor can be cancelled only by 
the addition of a specific homogeneous solution to (|2.8|) . States which are regular on both the past 
and future horizon have p^ + q'^ = 1, pq = Q and zero flux. These conditions give the Bunch-Davies 
state with the non- vanishing de Sitter invariant stress tensor, T^* = T^l = NH'^/24tt. The analog 
of the Boulware state with p = q = diverges as / — > exactly as in ()2.36|) with Th = H/2tt. 

In all the two dimensional examples considered, the parameter q may be viewed as an effective 
topological charge carried by the auxiliary field of the trace anomaly. The classical solutions of this 
order parameter field generally depend on global coordinate invariant quantities, such as \n{—K°'Ka) 
which diverges on a horizon where a timelike Killing field i^" becomes null, notwithstanding the 
finiteness of the local curvature there. Thus the scalar auxiliary field provides a coordinate invariant 
characterization of the local divergences of the stress-energy on the horizon, in macroscopic quantum 
states defined globally in the spacetime. 

III. EFFECTIVE ACTION AND STRESS TENSOR IN FOUR DIMENSIONS 

The low energy effective action for gravity in four dimensions contains first of all, the local 
terms constructed from the Riemann curvature tensor and its derivatives and contractions up to 
and including dimension four. This includes the usual Einstein-Hilbert action of general relativity, 

SEH[g] = -^ I d^x,/^{R - 2A) (3.1) 

as well as the spacetime integrals of the fourth order curvature invariants, 

silM = 11 V^9{aCataiC^'"^' + m')d'x, (3.2) 
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with arbitrary dimensionless coefficients a and /3. There are two additional fourth order invariants, 
namely E =*Rabcd*R"'^'^'^ and D-R, which could be added to H3.2|) as well, but as they are total 
derivatives yielding only a surface term and no local variation, we omit them. All the possible local 
terms in the effective action may be written as the sum, 

Siocai [9] = Y^ d'x ^^g{R- 2A) + 5^1 + Y. ^S ■ (3-3) 

n=3 

with the terms in the sum with n > 3 composed of integrals of local curvature invariants with 
dimension 2n > 6, and suppressed by M~^"+^ at energies E <C M. Here M is the ultraviolet 
cutoff scale of the low energy effective theory which we may take to be of order Mpi. The higher 
derivative terms with n > 3 are irrelevant operators in the infrared, scaling with negative powers of 
6°^° under global rescalings of the metric, as in (|2.2|) . and may be neglected at macroscopic distance 
scales. On the other hand the two terms in the Einstein-Hilbert action n = 0, 1 scale as e^"° or e^"'^ 
respectively, and are clearly relevant in the infrared. The fourth order terms in 1)3. 2() are neutral 
under such global rescalings. 

The exact quantum effective action also contains non-local terms in general. All possible terms 
in the effective action (local or not) can be classified according to how they respond to global Weyl 
rescalings of the metric. If the non-local terms are non-invariant under global rescalings, then they 
scale either positively or negatively under gab — > e'^'^^Oab- If "i^^ is some fixed length scale associated 
with the non-locality, arising for example by the integrating out of fluctuations of flelds with mass 
m, then at much larger macroscopic distances (mL ^ 1) the non-local terms in the effective action 
become approximately local. The terms which scale with positive powers of 6°^° are constrained 
by general covariance to be of the same form as the n = 0, 1 Einstein-Hilbert terms in SiocaU 
(|3.1I) . Terms which scale negatively with e^° become negligibly small as mL ^ 1 and are infrared 
irrelevant at macroscopic distances. This is the expected decoupling of short distance degrees of 
freedom in an effective field theory description, which are verified in detailed calculations of loops 
in massive field theories in curved space. The only possibility for contributions to the effective 
field theory of gravity at macroscopic distances, which are not contained in the local expansion of 
(|3.3j) arise from fluctuations not associated with any flnite length scale, i.e. m = 0. These are the 
non-local contributions to the low energy EFT which include those associated with the anomaly. 

Classical fields satisfying wave equations with zero mass, which are invariant under conformal 
transformations of the spacetime metric, Qab -^ s^'^dab have stress tensors with zero classical trace, 
T^"" = 0. Because the corresponding quantum theory requires a UV regulator, classical conformal 
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invariance cannot be maintained at the quantum level. The trace of the stress tensor is generally 
non-zero when H ^ 0, and any UV regulator which preserves the covariant conservation of T^ (a 
necessary requirement of any theory respecting general coordinate invariance) yields an expectation 
value of the quantum stress tensor with the non-zero trace 



{T^'') = bF + b'(E-^nR] +b"nR + ^PiHi, 



(3.4) 



in a general four dimensional curved spacetime. This is the four dimensional analog of H2.1() in two 
dimensions. In Eq. 1)3. 4|) we employ the notation, 



E =*Rabcd *R"' '^ — Rabcd.R"' '^ — '^RabR"' + R ; and (3.5a) 

with Rabcd the Riemann curvature tensor, Rabcd = ^abefR^ cd/'^ ^^^ dual, and Cabcd the Weyl con- 
formal tensor. The coefficients 6, b' and b" are dimensionless parameters proportional to h. These 
terms appear in the trace of stress tensor of any massless field, provided only that it is regularized 
and defined in such a way as to preserve general coordinate invariance and the covariant conserva- 
tion law, V(j(T'^) = 0, which we assume if the theory is to preserve general covariance. Additional 
terms denoted by the sum ^^ (iiHi in ()3.4() may also appear in the general form of the trace 
anomaly, if the massless field in question couples to additional long range gauge fields. Thus in the 
case of massless fermions coupled to a background gauge field, the invariant H =tr{FabF ) will 
appear in 1)3.4(1 with a coefficient /? determined by the beta function of the relevant gauge coupling. 
We shall concentrate on the universal gravitational terms in the anomaly in most of what follows, 
returning to consider the effect of the '^- PiHi terms at the end of this section. 

The form of (|3.4() and coefficients b and b' do not depend on the state in which the expectation 
value of the stress tensor is computed. Instead they determined only by the number of massless 
fields |g|, 

b = Y^^^ (Ns + QNf + 12Nv) , (3.6a) 

b' = - ^qqUttV (^^ + 1^^^ + 62^^) ' (^•6'") 

with {Ns, Np, Ny) the number of fields of spin (0, 2 , 1) respectively and we have taken h = 1. 

Three local fourth order curvature invariants E,F and OR appear in the trace of the stress 
tensor, but only the first two (6 and b') terms of (|3.4p cannot be derived from a local effective action 
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of the metric alone. If these terms could be derived from a local gravitational action we could simply 
make the necessary finite redefinition of the corresponding local counterterms to remove them from 
the trace, in which case the trace would no longer be non-zero or anomalous. This redefinition of a 
local counterterm (namely, the R^ term in the effective action) is possible only with respect to the 
third b" coefficient in ()3.4() . which is therefore regularization dependent and not part of the true 
anomaly. Only the non-local effective action corresponding to the b and b' terms in (|3.4|) lead to 
the possibility of effects that extend over arbitrarily large, macroscopic distances, unsuppressed by 
any ultraviolet cutoff scale. The distinction of the two kinds of terms in the effective action (local 

n 

or not) is emphasized in the cohomological approach to the trace anomaly |10||. 

The number of massless fields of each spin, Ns, Np, Ny is a property of the low energy effective 
description of matter, having no direct connection with physics at the ultrashort Planck scale. 
Indeed massless fields fluctuate at all distance scales and do not decouple in the far infrared. As in 
the case of the chiral anomaly with massless quarks, the b and b' terms in the trace anomaly were 
calculated originally by techniques usually associated with UV regularization (such as dimensional 
regularization, point splitting or heat kernel techniques) .9]. However just as in the case of the 
chiral anomaly in QCD, 1)3. 4() and 1)3. 6(1 can have significant effects in the far infrared as well. 

To find the WZ effective action corresponding to the b and b' terms in (|3.4)) . introduce as in two 
dimensions the conformal parameterization (|2.14|) . and compute 

V^F = V^F (3.7a) 

V^ U-^nR] = V^ (E-'^nll\+iVS^4'T, (3.7b) 

whose a dependence is no more than linear. The fourth order differential operator appearing in 



this expression is 1£ 
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A4 = n' + 2R^'VaVb - ^Rn + ^(V°i?)V, , (3.^ 



which is the unique fourth order scalar operator that is conformally covariant, viz. 

V^A4 = V^A^ , (3.9) 

for arbitrary smooth cr(x) in four (and only four) dimensions. Thus multiplying 1)3. 4|) by ^/—g and 
recognizing that the result is the a variation of an effective action Tyyz, we find immediately that 
up to terms independent of a, this effective action is 



WZ 



[g;a] = b (fxyf^Fa + b' I (fx.f^ | (e--Ur\ o- + 2o-A4Cj1 . (3.10) 
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This Wess-Zumino action is a one-form representative of the non-trivial cohomology of the local 
Weyl group and satisfies p.5|) just as its two dimensional analog T\yz does. Explicitly, by solving 
1)3 -Tbl) formally for a and substituting the result in p.lOf) we obtain 

'^Wz[g; a] = A^o Sanom[g] = Sanom[g = (?'^g\ - SanomM, (3-11) 

with 



Sanom[9] = ^ fd^'x./^ fd'^x'^/^fE-^aR] A^\x,x') 2bF + b' (e - ^U r) 



(3.12) 



and A^ (x, x') denoting the Green's function inverse of the fourth order differential operator defined 
by (|3.8)) . From the foregoing construction it is clear that if there are additional Weyl invariant 
terms in the anomaly (|3.4|) they should be included in the Sanom by making the replacement 
6F — > 6F + Yl,i PiHi in the last square bracket of Eq. ()3.12|) . 

Because of ()3.11() the anomalous effective action is clearly fixed only up to an arbitrary (local 
or non-local) Weyl invariant functional Sinv[g] of the metric, obeying 

A^ o Sinv[g] = Sinvie^^'g] - Sinvig] = . (3.13) 

Hence the general non-local form of the exact effective action of gravity must be of the form, 
Sanom + Sinv, and bchavior under the local Weyl group enables us to classify all terms in the exact 
effective action into three parts, viz. 



^ exact 



'-'local ~r ^anom ~r Jinv j V*^ •■'-'*/ 



given by (|3.3j) . (|3.12j) and (|3.13j) respectively. As in two dimensions, the anomalous term (|3.1()|) or 
H3.12() scales linearly with ctq, i.e. logarithmically with distance, and is therefore marginally rele- 
vant in the infrared with respect to the usual classical flat space fixed point [lO|. The fiuctuations 
generated by Sanom also define a non-perturbative Gaussian infrared fixed poin t, with conformal 
field theory anomalous dimensions analogous to the two dimensional case [2j|. This is possible 
only because new low energy degrees of freedom are contained in Sanom which can fiuctuate inde- 
pendently of the local metric degrees of freedom in Seh- Thus the effective action of the anomaly 
Sanom should be retained in the EFT of low energy gravity, which is specified then by the first 
two strictly relevant local terms of the classical Einstein-Hilbert action 1)3. H) . and the logarithmic, 
marginally relevant Sanom, i-e. 

Seffig] = SEH[g] + Sanom[g] (3.15) 
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contains all the infrared relevant terms in low energy gravity for E <^ Mpi. 

The low energy (Wilson) effective action (|3.15j) . in which infrared irrelevant terms are system- 
atically neglected in the renormalization group program of critical phenomena is to be contrasted 
with the exact (field theoretic) effective action of (|3.14jl . in which the effects of all scales are in- 
cluded in principle, at least in the approximation in which spacetime can be treated as a continuous 
manifold. Ordinarily, i.e. absent anomalies, the Wilson effective action should contain only local 
infrared relevant terms consistent with symmetry i23] • However, like the anomalous effective action 
generated by the chiral anomaly in QCD, the non-local Sanom niust be included in the low energy 
EFT to account for the anomalous Ward identities, even in the zero momentum limit; and indeed 
logarithmic scaling with distance indicates that Sanom is an infrared relevant term. Also even if no 
massless matter fields are assumed, the quantum fluctuations of the metric itself will generate a 



term of the same form as Sn 



0. 



Since all the infrared relevant terms are now included in Sg//, 
it contains all the low energy degrees of freedom of gravity, including also those massless scalar 
degrees of freedom in the local auxiliary field description which we exhibit below, and which are 
those responsible for the non-perturbative infrared fixed point with vanishing A found in 22]. In 
contrast, the additional terms in Sexact, being either marginally or strictly irrelevant in the infrared 
limit, contain non-universal and state dependent contributions which are not required for a consis- 
tent EFT expansion of the gravitational effective action in powers of E/Mpi. Since anomalies are 
the only way known that the decoupling hypothesis can fail, the neglected terms in Smv just as 
those in Siocai for n > 2 are also not expected to require the addition of any new degrees of freedom 
of the gravitational field which are relevant at low energies or macroscopic distances. 

The variations of the two {b and b') terms in Sanom give rise to two new conserved tensors in 
the equations of low energy gravity. To find their explicit form and exhibit the new scalar degrees 
of freedom they contain, it is convenient as in the two dimensional case to rewrite the non-local 
action Sanom in local form by introducing auxiliary fields. Two scalar auxiliary fields satisfying 

Aiip = UE-'^nR] , (3.16a) 

A4^|J = ^F, (3.16b) 

respectively may be introduced, corresponding to the two non-trivial cocycles of the b and b' terms 
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in the anomaly [li|. It is then easy to see that 

Sanom[g;V,'4'] = "^ / d^X y/^ < -(f/\iip+ {E- I^^R) f \ 

+ \ j d^x./^ [-2^l^^tl; + F^+(^E-'^nR\tl\ (3.17) 

is the desired local form of the anomalous action ()3.12|) ,,24,]. Indeed the variation of H3.17() with 
respect to the auxiliary fields (/? and ijj yields their Eqs. of motion H3.16() . which may be solved for 
if and ip by introducing the Green's function Aj {x,x'). Substituting this formal solution for the 
auxiliary fields into (|3.17|) returns (|3.12|) . The local auxiliary field form (|3.17|) is the most useful and 
explicitly contains two new scalar fields satisfying the massless fourth order wave equations H3.16() 
with fourth order curvature invariants as sources. The freedom to add homogeneous solutions to ip 
and ip corresponds to the freedom to define different Green's functions inverses A J (x, x') in (|3.12() . 
The auxilliary scalar fields are new local massless degrees of freedom of four dimesnional gravity, 
not contained in the Einstein-Hilbert action. 

As in the two dimensional case, the four dimensional anomalous action has an invariance, in 
this case "0 — *■ V' + ^o j with the corresponding Noether current given by: 

jab ^^ ab ' " -"- ^" ^ • 



r = V''U^ + 2[W^- -g"" Vbf - -n" + -VR , (3.18) 

where r^" is the topological current whose divergence, 

V„0» = E = RabcdR"'"'^ - ^RabR""^ + R^ (3.19) 

is the Euler-Gauss-Bonnet integrand in four dimensions. The current J" is conserved, 

VaJ°' = ^iV-- + \uR = Q (3.20) 

by (|3.16a|) . and the conserved charge may be taken to be 



o / J"dSa (3.21) 

167r2 I " "■ ^ 



with S a spacelike Cauchy surface on which initial data are specified or 

lav 



16vr2 '^ ■^'^ "^^^ (^-^^^ 



with dV the boundary of a Euclidean four volume V. The normalization of the Noether charge is 
chosen so that q = 1 corresponds to a delta function source of unit strength in the four dimensional 
Euler number, 

X = ^ f d^xV^E. (3.23) 
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The constant shift m the other auxihary field, i.e. ip ^ ip -\- ipo generates a global rescaling of 
the metric, and is not a symmetry of the action (|3.17|) unless E — 2\I\R/3 = F = 0. It is worth 
remarking that this enhanced conformal symmetry applies to solutions of the vacuum Einstein Eqs. 
if and only if the cosmological term A = 0. That the effective action of the anomaly may provide 
a mechanism to relax the cosmological term to zero was first proposed in ref. [2^. 

By using the definition of A4 and integrating by parts, we may express the anomalous action 
also in the form, 

Janom = b ^anom ' "^anom : I'^-'^^J 

with 

sZL^ J d^x,/^ |_(n^)(nV;) + 2(i?«^-|5'^'')(v,^)(v,V') 

+ ^Fip + ^(^E-'^nR^i;^ (3.25) 

We note that the S^^' and S^^' terms are very similar. Each is composed of terms at most quadratic 
and linear in the auxiliary fields ip and ip. Defining the general quadratic action, 

A[g■,^,^l;] ^ j d'^x./^g |- {U^) (D^) + 2 {r^^ - ^g^'^ (Va(/.)(V5^)| (3.26) 

and the two linear actions, 

B[g-^] = ]^j<fx./^{E-'^UR\^; (3.27) 

C[g;ip] = ^J d^xV^F^, (3.28) 

as functionals of the metric and the two auxiliary scalar fields (p and ip, we can express (|3.25p in 
the form, 

Sinlm = \ ^[5; V', V] + B[g; ^] ; (3.29a) 

Si^L = A[g- (^, V] + B[g- ^] + C[g; ^] . (3.29b) 

This implies that the covariantly conserved stress tensors derived from the E and F terms in the 
effective action, namely, 

Eab = — ^^ . ab (3-30a) 



-9 Sg" 

uoan 



2 SS'-^^ 
Fab = - ^- . ,, • (3.30b) 
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can be obtained from one fundamental tensor quadratic in the scalar auxiliary fields, 

+2 V^ [(Vc99)(VaVbV') + (VeV)(VaVbVP)] - ^ VaVfe [(VcV^)(V^^)] 

+^ i?ab (V,^)(V^V') - 4i2% [(V,)(^)(VeV) + (V,)V)(Vc(/9)] + ^ i?(V(,^)(Vfe)V') 
+^<7a6 {-3(n(^)(nV) + □ [(V,(^)(V^V')] +2 (Si?^'^ - i?5''') (Vc^)(VdV')} , (3.31) 

which is traceless, 

A:\g- ^,^]=^, (3.32) 

for arbitrary (p and ^, and two fundamental tensors linear in the scalar auxiliary fields, 



B,,[g;ip] = - A^" ;Y = -^V.V,n^-4C,Vv,V,9^ 



2_ ,5^[g;(/.] _ 2 

5 ^g"^ 3 

8 4 2 

+ ^<7ab {2nV + 6i?"'^VcVd(^-4i2n(/J + (V^i?)VcV7} , (3.33) 



Ca,[g;^] ^ -^^%^ = -4VeV, (qA)V) - 2 C^^ii'^'^y. . (3.34) 



and 



The latter of these is also traceless, whereas the Bab tensor has a non-zero trace whose value depends 
on which auxiliary field is substituted, i.e. 

Ba^ig; ^] = 2A4V9 = E - -UR (3.35a) 

5„"[5;V] = 2A^i; = F. (3.35b) 

From (HT^ . ^2B, (E2ni), (ESU, (pr:^ and (|T!ni) we have 

rpianom) _ 2 oSgnom _ L/p , t 77 /q qc^ 



-9 ^9' 

with 

Eab = -^Aabig; if, ^] + Bab[g; ^p] and (3.37a) 

Fab = Aab [g; ^, i^] + Bab [s; "0] + Cab [g] ^] ■ (3.37b) 
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Explicitly, 



Eab = -2(V(,(^)(Vfe)n(^) + 2V^[(V,v^)(V„Vfe(^)]-|v,Vfe[(V,(p)(VV)] 



3 

+ ^ Rab (VcV?)(VV) - 4i?^(JVfe)VP)(Ve(/^) + ^R (Va(/j)(Vb(^) 



2 „ _ . __ . ^, _ . _ . 2 



+ 6 



i^ab {-3(nv^)2 + D [(Ve^)(VV)] + 2 (Si?'^'' - i?/"^) (Ve(^)(V,(^)} 



and 



- VaVbD^ - 4C„V VcVdV? - 4i?[^Vfe)VcVJ +-RabU^ + -RVa^b^ 

-\{y(aR)'^b)^ + \9ab[2U^V + QR"^Vc'^d^-^RU^ + {y'R)yc^] , (3.38) 



Fab = -2 (V(,v?)(V;,) DV') - 2 (V(,^)(Vb) n(/^) + 2 V^ [(VeV?)(VaVW^) + (VeV')(V„V;,vp)] 
-^ V,Vb [(V,^)(V^V')] + ^ i?afe (V,(/?)(V^V) - 4i?% [(Vfe)(^)(V,V) + (V,)^)(Ve(/^)] 

+^i2(V(,</.)(Vb)V) + ^5a6 { - 3(n(/?)(n^) + D [(Ve(^)(v^v)] 

+2 {SR'^'^ - Rg^'^) {V,^){V,^l,)} - 4VeV, (C(,%V) - 2C7,V^^V 
-^ VaVbDV' - 4^'=,'^ V.VdV - 4i?^JVb)VeV) + \ Rab □V' + ^RVaVbi^ 



2 
3 



;V(,i?)Vb)V + ^5afe {2nV + 6i?^'^VeVdV-4i?n^ + (V^i?)(VeV')} ■ (3.39) 



Each of these two tensors are individually conserved and they have the local traces, 

E^" = 2A4Lp = E- -DR, (3.40a) 

o 

F/ = 2A,^ = F = CabcdC"'"" , (3.40b) 

corresponding to the two terms respectively in the trace anomaly in four dimensions (with /3j = 0). 
Results for the Eab tensor in terms of the decomposition (|2.2D , and the Fab tensor in a slightly 



different notation were presented previously in refs. 25| and [l0|, and partial results for Eab appear 
also in |2J]. 

It is straightforward to generalize the above development to include the effects of any additional 
Weyl invariant terms in the trace anomaly, transforming under ()2.2() by 

^f^Hi = ^H,. (3.41) 

By making the replacement F ^ F + Ylit^i^i/^ i^ (|3.16a|) . p.l7() . and (|3.25j) . we observe that 
C[g; ip] should undergo the replacement, 

C[g- V9] ^ C[g; ^] +Y^^^ J d'x V^ H,^ , (3.42) 
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with the corresponding additional terms in Cab and Fab in (jSSJ, ()3.37b|) . (|mffl|) . and (|3.40b|) to 
take account of any additional contributions to the trace anomaly. Hence no additional scalar 
auxiliary fields beyond ip and ^ are required in the general case. 

As in two dimensions, Sanom in four dimensions contains massless scalar degrees of freedom 
which remain relevant at low energies and macroscopic distances, and must be retained in addition 
to the usual EFT local derivative expansion in order to take account of the non-decoupling of 
massless fields. In contrast, the Weyl invariant terms in ()3.2() or Sinv do not contain additional 
degrees of freedom of the gravitational field which are relevant at energies far below the Planck 
scale Mpi^ and are either marginally or strictly irrelevant in the infrared. 

In a given fixed spacetime background the solutions of the fourth order differential equations 
1)3. 16(1 for the scalar auxiliary fields (/? and tj: may be found, and the results substituted into the 
stress tensors Eab and Fab- The freedom to add homogeneous solutions of (|3.16|) to any given inho- 
mogeneous solution corresponds to the freedom to change the Weyl invariant part of the effective 
action and corresponding traceless parts of the stress tensors, without altering their local traces 
H3.4U() . The traceless terms in the stress tensor thus depend on the state of the underlying quantum 
field(s) which can be fixed by specifying boundary conditions on the auxiliary fields. When the 
spacetime is conformally flat or approximately so, the Weyl invariant action Smv can be taken to 
vanish, as it does in the conformally related flat spacetime, and the low energy effective action 
1)3. 15(1 should become a good approximation to the exact effective action 1)3.14(1 . In order to test 
this hypothesis we evaluate the corresponding anomalous stress tensor (|3.36|) in conformally flat 
spacetimes for a number of different macroscopic states. 

IV. STRESS-ENERGY TENSOR IN CONFORMALLY FLAT SPACETIMES 

A. Flat Spacetime 

In flat Minkowski spacetime, E = F = n-R = 0, and the auxiliary fields satisfy the homogeneous 
wave eqs.. 



A, 



(^= nV = o= D^V'- (4.1) 

flat 



For the trivial solutions (/9 = -0 = (or an arbitrary spacetime constant), the stress tensor van- 
ishes identically. This corresponds to the ordinary Lorentz invariant vacuum state, with vanishing 
renormalized {TJ'). Hence flat spacetime will continue to satisfy the Einstein equations modified 
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by the trace anomaly terms (with vanishing cosmological term). 

In flat spacetime the wave Eq. (|4.H) has plane wave solutions obeying the massless dispersion 
relation, uji^ = ±|k|. Since the wave operator here, D^ is fourth order, there are two sets of 
massless scalar modes for each auxiliary field, (p and ^. These are new local scalar degrees of 
freedom for the gravitational field not present in classical general relativity, and give rise to scalar 
gravitational waves. These scalar gravitational waves of the augmented effective theory have as 
their sources the fourth order curvature invariants, E,F and \3R which are quite small for most 
astrophysical sources, excepting perhaps those of cosmological origin. Hence in the flat space limit 
the scalar degrees of freedom decouple completely. Even in the presence of matter, they couple 
only indirectly through the effects of spacetime curvature, and only very weakly at that. Localized 
sources act as only very weak generators of scalar gravitational radiation in the (p and V fields, 
and primordial scalar gravitational radiation, even if present, would be difficult to detect, except 
indirectly through the gravitational effects of the stress-energy H3.36|) . The effects on the production 
of scalar gravitational waves from the coupling of ip to gauge fields through the Hi terms in the 
trace (|3.4() is also possible, and merits a separate investigation. 

For non-vanishing ip the tensor Eat in flat space becomes 

Eab = -2 (V(„^)(Vb) DV?) + 2{n^){VaVb^) + -(Vc(/j)(V^VaVfe(^) - - (VaVeV?)(V6VV) 

flat o o 

+7:9ab {-3(n^)2 + n(VcV9VV)}-^VaVbn(/p, (4.2) 

3 

in general curvilinear coordinates. We note that although the energy density for the scalar auxiliary 
fields is not positive definite, this does not lead to any instability at macroscopic distances in flat 
space, because of the decoupling of the auxiliary fields in the flat space limit. In ref. [251 the physical 
state space of the conformal scalar fluctuations of Sanom (proportional to 93), once decoupled from 
metric fluctuations, was shown to have positive norm, and be free of any tachyon or ghost modes. 
When Sanom is considered together with the classical Einstein-Hilbert terms, because of the higher 
derivatives in the auxiliary field stress tensor, instability in flat space appears only at ultrashort 
Planck scale wavelengths 26|, i.e. at the order of the cutoff scale and outside the range of validity 
of the low energy EFT approach. At the Planck scale the local expansion of the effective action 
in powers of E"^ /Mp^ clearly breaks down. For E <C Mpi where the low energy effective action 
1)3. 15() can be applied, the auxiliary fields are coupled very weakly and indirectly to matter and 
propagate as very nearly non-interacting free field scalar excitations. This shows that empty flat 
space is stable to the effects of Sanom on macroscopic distance scales, and more subtle signatures 
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of the presence of Sanom in nearly flat spacetimes must be sought. 

As a simple example of how the stress tensor (j4.2j) may be used to evaluate macroscopic semi- 
classical effects of quantum fields in flat space, consider the special solution to (|4.1() . 

for some constant a with dimensions of length and ci dimensionless. With a similar ansatz for 
■0 = C2Z^ jlc? the total stress tensor of the anomaly (|.S.36|) becomes 



7^ar"^^ = :7^diag (-1,1, 1,-3), (4.4) 



with C = cfb' + 2bc2- This is the form of the stress tensor of the Casimir effect in region between 
two infinite parallel plates a distance a apart in the z direction. The constants ci, C2 and C depend 
on the boundary conditions which the quantum field(s) obey on the plates, and parameterize the 
discontinuities in the components of the stress tensor there, but the generic form of the stress tensor 
appropriate for the geometry of two parallel conducting plates follows from the simple ansatz ()4.3() 
for the classical auxiliary fields together with (|3.36)) . Thus the low energy effective action of the 
trace anomaly in curved space correctly determines the form of the traceless components of the 
macroscopic quantum stress tensor for the parallel plate geometry in flat space. 
If instead of (|4.3|) we choose 

^='JTh' (4.5) 

with Tp = C2T^t^/2 the stress tensor ()3.36() becomes 

ri~) = C^diag(-3,l,l,l), (4.6) 



with C = —cfb' + 2bc2 which is the form of the stress tensor for massless radiation at temperature 
T, which is again traceless in fiat space. The same form for the anomalous stress tensor is also 
obtained by considering the general static spherically symmetric ansatz, ip = ip{r). Then 

f = -^— + Co + cir + C2r^ , (4.7) 

r 

with c_i = to ensure a regular stress tensor at the origin. We note that although a thermal state 
is a mixed state for the radiation field, it is a macroscopic pure coherent state for the gravitational 
auxiliary potentials at the semi-classical level, where fluctuations around the mean stress tensor 
(r.,)=Ti™^ are neglected. 
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For the four dimensional Rindler metric, 

ds"^ = -p^drj^ + dp^ + dy'^ + dz'^ (4.8) 

and (p = (p{p), A^ip = O^ip = has the non-trivial solutions p^ In p,p^ and In p. The first two 
solutions lead to stress tensors which are non-vanishing at infinity similar to the thermal state, 
while the Inp solution for ip and ip gives 

r„'' = -^diag(-3, 1,1,1). (4.9) 

Adding a term linear in rj to the solution for <p{p) gives a stress tensor of the same form as (|4.9|) 
with a different constant C. In either case the form of the stress tensor ()4.9|) is that of the Fulling- 
Rindler state in four dimensions '27J. Because of 1)2.24(1 . the solution gives rise to a singularity in 
Dip at p = proportional to 6{p)/p, and therefore is a solution only in the distributional sense, as in 
the two dimensional case. The higher order derivatives in the Noether current in four dimensions 
produce higher order singularities at p = 0, and the conformally transformed metric e~^ds'^ is 
singular at p = for any ci / 0. 

B. de Sitter Spacetime 

In conformally flat spacetimes with Qab = e^'^^abj one can choose p = 2a and ip = to obtain 
the stress tensor of the state conformally transformed from the Minkowski vacuum. In this state 
Fgh vanishes, and p can be eliminated completely in terms of the Ricci tensor with the result 

lij, yy, 

Eab = —2 Wab — - 'Hab 

y 

= -Va^bR + 2R,'Rbc - -^RRab + gab ( "g ^^ " RcdR""" + g^' j • (4-10) 

Thus all non-local dependence on boundary conditions of the auxiliary fields y? and ^l^ vanishes in 
conformally flat spacetimes for the state conformally mapped from the Minkowski vacuum. In the 
special case of maximally 0(4, 1) symmetric de Sitter spacetime, Rab = "iH^dab with R = 12H^ a 
constant, the tensor ^^'Hab vanishes identically and 



E, 



ab 



= -2 ^'^H^b 
dS 



= &H^9ab . (4.11) 

dS 
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Hence we obtain immediately the expectation value of the stress tensor of a massless conformal 
field of any spin in the Bunch-Davies state in de Sitter spacetime, 



T, 



ab 



= 6b' H"" gab = 7T Qab {Ns + lliVf + 62 A^^) , (4.12) 

BDAS 9607r2 ^""^^ "^ / ^ ^^ ' ^ ' 



which is determined completely by the trace anomaly. 

It is also possible to consider states which are not maximally 0(4, 1) symmetric. For example, 
if de Sitter spacetime is expressed in the spatially flat coordinates, 



ds^ 



= -dr^ + e^^^dx^ (4.13) 

dS 



and if = C/3(t), we obtain from (|3.36() the stress-energy in spatially homogeneous, isotropic states. 
Since in this maximally symmetric spacetime the A4 operator factorizes, 

^={U-2H^){U^) = 12H\ (4.14) 

da 

it is straightforward to show that the general solution to this equation with (p = (p{t) is 

ip{T) = 2Ht + Co + c_ie-^" + C-ae-^^" + C-se'^^" . (4.15) 

The first (inhomogeneous) term can be understood geometrically from the fact that the conformal 
transformation, 



e-'^BD ds'^ 



e-2^^ds2 



dS 

1„-Ht 



= -dr]^ + dx^ (4.16) 

dS 



brings de Sitter spacetime to flat spacetime, with r/ = —He~ the conformal time. This value 
of 99 = 'Pbd{'t) = 2Ht is exactly the one that gives the Bunch-Davies stress-energy (|4.12|) when 
substituted into 1)3. 38() . When the full solution for c/?(r) of 1)4. 15() is substituted into ()3.38|1 we 
obtain additional terms in the stress tensor which are not de Sitter invariant, but which fall off 
at large r, as e~ . The stress tensor of this time behavior is traceless and corresponds to the 
redshift of massless modes with the equation of state, p = p/3. This is a special case of the detailed 
analysis of the stress tensor expectation value of a field of arbitrary mass in homogeneous, isotropic 
states in de Sitter spacetime of ref. i29|. It is straightforward to generalize these considerations 
to spatially homogeneous and isotropic states in arbitrary conformally flat cosmological spacetimes 
with Robertson- Walker scale factor a{r]). Because of the reduction of the auxiliary field stress tensor 
Eab to the local geometric form ()4.1U() in conformally flat spacetimes, it can have no significant effects 
in states with Robertson- Walker symmetries when the curvature is much less than the Planck scale. 
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States of lower symmetry in de Sitter spacetime may be considered by taking static coordinates, 



dr 



ds' = -f{r)dt' + T7T + r'dn' , 
fir) 



(4.17) 



where 



f{r] 



dS 



1 - ijV^ , 



(4.18) 



If we insert the ansatz ip = (p{r) in (|4.14|) . the general 0(3) spherically symmetric solution regular 
at the origin is easily found: 

l-Hr" 



V{r] 



ln(l-^V)+co + |ln(l^U^^..^ln 



(4.19) 



dS ^ ' '^ 2 \L + HrJ 2Hr \l + Hr ^ 

A possible homogeneous solution proportional to 1/r has been discarded, since it is singular at 
the origin. An arbitrary linear time dependence 2Hpt could also be added to '^{r), i.e. (p{r) — > 
(p{r,t) = ip{r) + 2Hpt. The particular solution. 



cp^^ (r, t) = In (1 - H^r^) + 2Ht = 2Ht . 



(4.20) 



is simply the previous solution for the Bunch-Davies state we found in the homogeneous flat coor- 
dinates ()4.13|) . The priveledged role of this fully 0(4, 1) invariant state may be understood in the 
static coordinates by recognizing that the line element 1)4. 17() may be put into the alternate forms, 



ds' 



f 



„*2 



dt' + dr*^ + —dn 



f [-dt^ + e^ {dx^ + sinh^ X dn"^)] , 
[-dp2 + p2 {dx^ + sinh^ x dn^)] 



i'f 



4 i-dr]^ + dx^) , 
e'^^^ids^)fiat. 



by means of the successive changes of variables, i = H ^, 

>"■ dr 



(4.21) 



Hr* =H 



m 



ilnfi±^ 
2 \1-Hr 



X 

t = -£lnp, 

r] = ip cosh X ) 

X = ip sinh X (sin 9 cos (j) , sin 9 sin (j) , cos 



(4.22) 
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Thus de Sitter spacetime in static coordinates is conformally related to flat spacetime by the 
conformal transformation e^^^ with ipsD given by H4.2U() as before. The benefit of carrying this out 
exphcitly in the static coordinates is that the steps in (|4.21() can be repeated essentially unchanged 
in the vicinity of any static Killing horizon, without reliance on the special higher symmetries of 
de Sitter spacetime. Then 

^h{r,t)=\nf{r)±jt (4.23) 

is the value of the auxiliary field which conformally transforms the local neighborhood of the 
horizon to flat spacetime. In de Sitter space this leading order conformal transformation mapping 
the vicinity of the horizon to flat space is exact, c/. (4.21). 

In the general 0(3) symmetric state in de Sitter spacetime the radial component of the Noether 
current, 

r = v^ny?-2F2v>--J7^ 

in the frame in which Q''' = SH'^r. Thus the parameter q coincides with the definition (|3.21|) and if 
non-zero gives rise to a delta function singularity for A4ip at the origin r = 0. It is also the strength 
of the puncture configurations (also called "spikes") found in J30| . 

In the general 0(3) symmetric state centered about the origin of the coordinates r = 0, the 
stress tensor is generally dependent on r. In fact, it generally diverges as the observer horizon 
r = H~^ is approached. From (|4.19|) we observe that 

/I - Hr\ 
ip{r) ^c^ln — ]+0{l-Hr) as Hr ^ 1 , (4.25) 

dS \ 2 J 

so that the integration constant c^ controls the singular behavior at the observer horizon r = H^^. 
As in the Schwarzschild case a non-zero c^ produces a higher order singular distribution for A493 
at the horizon. 

The second auxiliary field ip satisfies the homogeneous Eq., A4'(/; = 0, which has the general 
spherically symmetric solution linear in t, 

. , , q' fl-Hr\ 2d„-q' f I - Hr\ 

Note that the constant d^ enters this expression differently than c^ enters the corresponding 
Eq. ()4.19|1 . due to the inhomogeneous term in 1)4. 14p . which is absent from the iIj equation. 
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Since the anomalous stress tensor is independent of cq and do, it depends on the six parame- 
ters {Cfj,dg,q,q' ,p,p') in the general stationary 0(3) invariant state. However, in order to have no 
singularity of the stress tensor at the origin we must choose q = q' = 0. With no sources or sinks 
at the origin, the choice q = q' = then satisfies the zero flux condition, 

Tj "" = -4H^ {bpq' + bp'q + b'pq) = (4.27) 

automatically, for any p and p'. For the other components of the stress tensor (|3.36|) . the conditions 
of finiteness on the observer horizon at r = H^^ are 

c^ = l, p = ±l, d„ = ±p' (4.28) 

which are satisfied by the Bunch-Davies state with p = 1 and d^ = p' = 0. On the other hand, the 
local vacuum state at the origin (i.e. the analog of the Boulware vacuum |3l| in de Sitter space) we 
should take p = p' = Q. Since this is inconsistent with 1)4. 28() . this state has a diverging stress-energy 



on the observer horizon. There is no analog of the Unruh state 33] in de Sitter spacetime, since by 
continuity a flux through the future or past observer horizon at r = H~^ would require a source or 
sink of flux at the origin r = 0, a possibility we have excluded by ()4.27|) above. 

V. NEAR HORIZON CONFORMAL SYMMETRY AND STRESS TENSOR 

The fact that the trace anomaly can account for even the tracefree parts of the stress tensor 
exactly in two dimensions and conformally flat spacetimes in higher dimensions is due to the fact 
that Sinv may be neglected in these cases. In these circumstances the low energy effective action 
1)3. 15(1 becomes a good approximation to the exact effective action (|3.14l) . However, as has been 
noted by several authors [331 , spacetimes with a static Killing horizon are locally conformally related 
to flat spacetime in the vicinity of the horizon. Because of this near horizon conformal symmetry we 
expect the anomalous action ()3.12|) to generate a stress tensor 1)3. 36(1 which is a good approximation 
to the quantum stress tensor near an arbitrary Killing horizon. 

To exhibit the conformal geometry of a static Killing horizon consider the general static spher- 
ically symmetric spacetime of the form, ()4.17|) where the function /(r) possesses a simple zero ar 
r = rj^, i.e., 

f{r)=±-(r-r+)+0{r-r+f . (5.1) 
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Defining r* in the usual way, 

dr 



f{r) 2 

where the elhpsis denotes regular terms as / — > 0, and 



±-ln/ + ..., (5.2) 



1 fff\ r* 
-In ^ * - 

2 V4r2 



X = -oln(74)=T- + ---, (5.3) 



we find that near r = r_|_, 



^2 /)2 

/(r) 4 

Hence we may repeat the steps (|4.22|) of the previous de Sitter case for an arbitrary spherically 
symmetric static Killing horizon, in the vicinity of the horizon. The length scale i is also related to 
the surface gravity of the horizon by £~^ = k. Thus we have the conformal mapping ()4.21|) of the 
near horizon metric (|4.17|) to the flat metric with 93 ~ fhi''', t) given by (|4.23|) . Since this conformal 
transformation is singular at r = r+ , the exact effective action (|3.14|) is dominated by its anomalous 
contribution Sanom, and the stress-energy is dominated by H3.36() in any state for which ip and the 
corresponding T^ diverges on the horizon. Hence the near horizon conformal symmetry (|4.21j) may 
be used to calculate the exact behavior of the divergent terms of the quantum stress-energy tensor 
from (|3.36|) in any spherically symmetric static spacetime with a Killing horizon. We illustrate this 
for the most important case of Schwarzschild spacetime. 

A. Schwarzschild Spacetime 

In the four dimensional Schwarzschild geometry, /(r) = 1 ~ r/2M — 1 -|- 0{r — 2M)'^. 

Hence r_(_ = 2M and i = AM. The conformal transformation that maps the near horizon geometry 
to flat spacetime is therefore 

V.;,(r,t)=ln/±^. (5.5) 

Since Schwarzschild spacetime is not globally conformally flat the conformal transformation of the 
near horizon geometry will receive subleading corrections to (|5.5|) . To find these subleading terms 
we require the general solution to 1)3. 16() . Noting that Schwarzschild spacetime is Ricci flat, we have 



F\^=E\^=RabcdR''''''' 



48 M^ 
V,0« = —^ . (5.6) 
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A particular solution of either of the inhomogeneous Eqs. (|l-{.16|) is given then by '^{r), with 

4M . / r \ 1 / AM' 



dip 
dr 



/ r \ i / 4M 

3r2/ '^ V2M; " 2M V^ ^ ~ 



(5.7) 



Notice that this particular inhomogeneous solution to (|3.1(ij) is regular as r — > 2M. The general 
solution of (|3.16|) for ip = ip(r) away from the singular points r = (0, 2M, oo) is easily found and 
may be expressed in the form |24| . 



dip 
dr 



dp 
dr 



+ 



2Mc, 



+ 



drr 



s rV 4M2r2/ 72^^, 

2 / r 2M\ ( 2M\ q 
77TT + 1 + Inl ]-^ 



■' 2M \2M r J 



6M \2M 
1 
~3M 



r J 

2Mc 



r J 6r 



4M , / ^ \ ^ 



H 



1 

r^ r{r-2M) ' 2M \2M 



Coo f r 2M 



(5.8) 



in terms of the three dimensionless constants of integration, c^j, Cqo, and q. This expression has 
the limits, 



dp 
dr 
dp 
dr 



2M 



+ 



2M "V2M 



C^rxn' 



+ 



2c, 



oo Q' Cqo 



-1 

2M 



1 

2M 



3Cn 



- ^-3' + 



■ g In 



(^\ 



+ 



2M 



18 



(9-2) 



r ^ 2M, (5.9a) 
+ ..., r ^ 00. (5.9b) 



4M2 4M r 3r2""""V2My 

Hence c^ controls the leading behavior as r approaches the horizon, while Coo controls the 
leading behavior as r — > 00, which is the same as in flat space. The leading behavior at the horizon 
is determined by the homogeneous solution to (|3.16|) . c^ In/ = c^ \n[—K"-Ka)'^ where K = dt is 
the timelike Killing field of the Schwarzschild geometry for r > 2M. To the general spherically 
symmetric static solution (|5.8j) we may add also a term linear in t, i.e. we replace p{r) by 



plr, t) = pir) -\ t . 



(5.10) 



Linear time dependence in the auxiliary fields is the only allowed time dependence that leads to a 
time-independent stress-energy. We see that the conformal transformation to flat space near the 
horizon 1)5. 5() corresponds to the particular choice c^j = zizp = 1, leaving the subdominant terms in 
(|5.9a|) parameterized by q and Cqo undetermined. 
The topological charge q may be identified from 

8M2 



VDp 



+ 



(5.11) 



r^ 2Mr2 ' 

so that q is given by ()3.22|) in the Euclidean signature Schwarzschild geometry in the frame in which 
0'" = — 16M2/r^ and O* = 0. A non-zero q gives A^p a delta function singularity at the origin. 
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Non-zero c^ or Coo produce higher order singular distributions in A^ip at the horizon or infinity 
respectively. We note that unlike in flat space there are logarithmic terms in (p{r) and no value of 
q which can eliminate the subleading logarithmic behavior of ip' at both r = 2M and r = oo. 

Since the equation for the second auxiliary field ■0 is identical to that for (p, its solution for 
■0 = ip{r,t) is of the same form as ()5.8|) and H5.1U() with four new integration constants, d^,doo, 
q' and p' replacing c^,Coo, q, and p in {p(r,t). Adding terms with any higher powers of t or more 
complicated t dependence produces a time dependent stress-energy tensor. Inspection of the stress 
tensor terms in (|3.36j) also shows that it does not depend on either a constant ipo or ^o but only the 
derivatives of both auxiliary fields in Ricci flat metrics such as Schwarzschild spacetime. For that 
reason we do not need an additional integration constant for either of the fourth order differential 
equations ()3.16|) . 

With the general spherically symmetric solution for ip{r,t) and %l){r,t), we can proceed to com- 
pute the stress-energy tensor in a stationary, spherically symmetric quantum state. For example 
the Boulware state may be characterized as that state which approaches the flat space vacuum as 
rapidly as possible as r ^ oo J3l|. In the flat space limit this means that the allowed r"^ and r 
behavior in the auxiliary fields (r and constant behavior in their first derivatives) must be set to 
zero. Inspection of the asymptotic form (|5.9b|) shows that this is achieved by requiring 

Coo = '^oo = and (5.12a) 

q = (i = 0. (Boulware) (5.12b) 

If we set p = p' = as well, in order to have a static anasatz for the Boulware state, then the 
remaining two constants c^ and d^ are free parameters of the auxiliary fields, and lead to a stress- 
energy which generically has the divergent behaviors, 

-2-1,2 11 r -2M . . 

s , s , In s, and ms, as s = — — >0, (5.13) 

on the horizon. It is not possible to cancel all four of these divergent behaviors with the two 
remaining free parameters. As in the two dimensional case, it is not possible to have auxiliary fields 
and stress-energies which both fall off at infinity and are regular on the horizon. This macroscopic 
effect of the trace anomaly is a result of the non-trivial global topology of the Schwarzschild metric, 
notwithstanding the smallness of local curvature invariants at the horizon. 

The linear time dependence 1)5. 1U() is consistent with the singular behaviors (|5.13|) of the stress- 
energy on the horizon. Since the stress-energy in the Boulware state diverges as s"^ in any case. 
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we can match this leading divergence in the anomalous stress tensor b y a djusting c^ and d^ 
appropriately. One then has a one parameter fit to the numerical data of 34] . The results of this 
fit of H3.36() with c^ and d^ treated as free paramters are illustrated in Figs. 1, for all three non-zero 
components of the stress tensor expectation value of a massless, conformally coupled scalar field in 
the Boulware state. The best fit values plotted were obtained with c^ = —^ and d^ = ft' 

For comparison purposes, we have plotted also the analytic approximation of Page, Brown, 



and Ottewill 



35 



36 



(dashed curves in Figs. 1). We observe that the two parameter fit with 



the anomalous stress tensor in terms of the auxiliary (p and ip fields is more accurate than the 



approximation of refs. 



35 



36l . l37| for the Boulware state. In the latter case the stress-tensor is 
approximated by making a special conformal transformation a = 2^^/ ~ 2^'^{~K"'^a) to the 

optical metric, obtained from the Schwarzschild line element 1)4. 17() by dividing by /(r) = 1 . 

This procedure was also motivated by the form of the trace anomaly, in that for ultrastatic metrics 
with the timelike Killing field if", the optical metric conformally related it has vanishing trace 
anomaly 36]. However, unlike the general form of the anomalous stress tensor in ()3.36() . which 



mvo. 



35, 



ves solving linear equations for two independent auxiliary fields, the approximation of refs. 
361 . l37| is applicable only to special cases, such as Schwarzschild geometry which are both 



ultrastatic and Ricci flat, and admit a particular conformal transformation to the static optical 
metric. In contrast the approach based on the auxiliary field anomalous effective action does 
not require any special properties of the background spacetime, and because of the near horizon 
conformal symmetry gives the general form of the divergences of the stress-energy on horizon. 
An important practical difference between these approximations and that considered here is that 
they yield simple algebraic approximations to the stress-energy, whereas ()3.36|) generally contains 
logarithmic terms as well. Evidently the general form of the effective stress-energy tensor due to 
the anomaly in terms of two auxiliary fields and the freedom to add homogeneous solutions to the 



linear Eqs. ()3.16() al 



than that of refs. 



35 



ows 



36, 



br a better approximation to the exact stress tensor in the Boulware state 



371 ] near the Schwarzschild horizon. 
The stress-energy diverges on the horizon in an entire family of states for generic values of 
the eight auxiliary field parameters {cg,q,Coo,P]dg,q',doo,p'), in addition to the Boulware state. 
Hence in the general allowed parameter space of spherically symmetric macroscopic states, horizon 
divergences of the stress-energy are quite generic, and not resticted to the Boulware state. On the 
other hand, the condition 1)2. 32() that the stress-energy om the horizon be finite gives four conditions 
on these eight parameters, in order to cancel the four possible divergences listed in 1)5. 13() . Referring 
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to the explicit form of the stress-energy components in static spherically symmetric spacetimes 
given in the Appendix, these four conditions can be satisfied in several different ways. The simplest 
possibility with the minimal number of conditions on the auxiliary field parameters is via 

{2b + b')cl+p{2bp' + b'p) = 0, (s-^), (5.14a) 

{b + b')c„ = bd^, (s-i), (5.14b) 

q = q' = 2 , (In s and Ins). (horizon finiteness) (5.14c) 

It is noteworthy that this cancellation requires a non-zero topological charge q = 2, much as in 
the two dimensional Schwarzschild case, and again unlike the Fulling-Rindler case in either two 
or four dimensions. The q = 2 condition may be understood from the fact that it cancels the 
logarithmic divergence of In / in the second form of (|5.8j) , leaving however a logarithmic singularity 
and non-trivial topological charge at both the origin and infinity. 

Also noteworthy is that the the auxiliary field configurations which lead to diverging stress 
tensors on the horizon have finite action (|3.24jl . This is clear from the fact that 

I^V' ^ TT79" 1^(7^77 ~ l) +co"'5*- as r^2M. (5.15) 



4M2 \2M 

Hence both {Otp)'^ and {\I\ip){\Ilip) are proportional to In (r/2M — 1), which is integrable with 
respect to the measure r'^dr as r — > 2M, and the integrals in H3.25() converge at the horizon for 
the general solution (|5.1()|) of the auxiliary field equations. This is not the case at infinity, unless 
the four constants Coo, doo,q and q' are all required to vanish. From (|5.8|) in this case the auxiliary 
fields fall off at large distances r ^ 2M similarly to the scalar potential of Newtonian gravity, or 
Brans-Dicke theory [33|. Because of the very different way they couple to spacetime curvature, the 
auxiliary fields of Sanom give only a very weak long-range interaction between massive bodies, that 
falls off very rapidly with distance (at least as fast as r ). Thus from consideration of the general 
form of the anomalous action and auxiliary fields in the vicinity of the Schwarzschild horizon and 
their falloff at infinity, we find no a priori justification for excluding generic states with stress-energy 
tensors that grow without bound as the horizon is approached. In such states the backreaction of 
the stress-energy on the geometry will be substantial in this region and lead to qualitatively new 
macroscopic effects on the horizon scale. 

If instead of requiring rapid falloff at infinity we require regularity of the stress-energy on the 
horizon, i.e. conditions 1)5. 14() . which is a four parameter restriction of the eight parameters in the 



auxiliary fields of the form (|5.8|) and 



Hawking-Israel and Unruh states 



3^, 




. These are the conditions appropriate to the Hartle- 
Because of the linear time dependence (|5.1U|) we also 
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obtain a no-zero flux (T^ '') in general. For the Hartle-Hawking state the vanishing of this flux gives 
a flfth condition, viz. 

b{qp' + q'p) + b'pq = , (zero flux) (5.16) 

on the parameters of the auxihary fields. Of the remaining three parameters two can be fit by the 
flnite values of (T^ *) on the horizon and at infinity respectively. This leaves one final parameter free 
to adjust for a best fit of the form of all three components of the stress-energy tensor at intermediate 
points. 

Alternately, the conformal transformation 1)5. 5|) could be used to flx the leading order behavior 
of the auxiliary field cp near the horizon, i.e. c^ = 1 = itp. Then 1)5. 16(1 fixes p' which with 
q = q' becomes redundant with the first of the conditions ()5.14l) . Hence only Cqo and doo are 
left undetermined, and no free parameters at all are left if the values of (Tj*) on the horizon and 
at infinity are fixed. Using the first method with one free parameter, the results for the (Tj *) 
component of the scalar, Dirac spin ^ and electromagnetic fields are shown in Figs. 2 (solid 
curves), for typical values of the parameters. The data points plotted in Figs. 2a and 2b are the 
numerical results of the direct quantum field theory calculations for spin and -^ of Teis. |40| and 
|4ll | respectively. In the spin 1 case the dashed curve represents the numerical results of ref. 42] 
in an accurate analytic polynomial fit to their data provided by the same authors. 

We note that the energy density in the Hartle-Hawking state computed from the general form 
of the anomaly is not in especially good agreement with the direct numerical calculations from the 
underlying quantum field theory. This is in accordance with the general discussion of the exact 
and low energy effective actions of the previous section. In states which have no divergences on the 
horizon there is no particular reason to neglect the Weyl invariant terms Sinv in the exact effective 
action, and the stress-energy tensor it produces would be expected to be comparable in magnitude 
to that from Sanom, and both contributions are negligibly small on macroscopic scales. However 
in states such as the previous Boulware example, the diverging behavior of the stress tensor near 
the horizon is captured accurately by the terms in (|3.36() arising from the anomaly, which have the 
same generically diverging behaviors as the quantum field theory expectation value {T^ )■ Even in 
states for which the expectation value remains regular on the horizon, the auxiliary field anomalous 
stress-energy can be matched accurately to the finite value of (T^'') on the horizon with suitable 
choice of integration constants in ()5.8|) . 

For comparison plotted also in Figs. 2a and 2b is the analytic approximation to (T^*) of Frolov 
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and Zel'nikov (FZ) [4j|. This latter approximation gives an improvement over that of |35l. ISfil. l37l| 

since it permits adjustment of the value of the stress-energy on the horizon to the numerical data, 
rather than fixing it to an incorrect value. Once that is done, the FZ approximation automatically 
gives the correct asymptotic value of the stress-energy in the Hartle-Hawking state at infinity, but 
contains no additional free parameters. For the spin ^ case the FZ approximation is a significant 
improvement of that of PBO which disagrees with the numerics even in the sign of (T^ *) on the hori- 
zon. Although it is a noticeably better fit to the numerical data than 1)3. 36(1 . the FZ approximation, 
like that of PBO, relies on the existence of a timelike Killing field and both are therefore highly 
specialized to certain spacetime backgrounds such as the Schwarzschild black hole. In contrast the 
stress-energy tensor 1)3. 36() arising from Sanom in the low energy EFT does not require the existence 
of a timelike Killing field, and can be computed in principle for arbitrary spacetimes. 

For the time asymmetric Unruh state the zero flux condition 1)5. 16(1 is replaced by the properly 
normalized finite fiux condition, 



W'>1 =-i^ P-17) 



where the luminosity, L is given by 



L=^ {bqP + bq'p + b'pq) (5.18) 

in terms of the auxiliary field parameters. Otherwise one can proceed as in the previous Hawking- 
Hartle state by fixing two of the three remaining auxiliary field parameters by the values of (T^ *) on 



the future horizon and at infinity which are known numerica 
fields in the Unruh state from the work of refs. |4J| and 45 



■ for the scalar and electromagnetic 
4a |. This leaves again one parameter 



free to fit. The results of that fit for the conformally coupled scalar and electromagnetic fields in 
the Unruh state are shown in Figs. 3 (solid lines). The dashed curves are accurate polynomial fits 
of [43] and 48 1 to the numerical data of [45| and |4a| for the spin and 1 fields respectively. There 
is at present no numerical calculation of the stress-energy tensor of the massless Dirac field in the 
Unruh state with which to compare. 

B. Other Spacetimes with Horizons 

The foregoing analysis of the anomalous action and stress tensor in the Schwarzschild back- 
ground may be extended to other spacetimes with metrics of the form (|4.17() . The electrically 
charged Reisner-Nordstrom black hole for which /(r) = 1 — 2M/r + e^/r^ is an example in which 
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Eqs. (|5.H) - (|5.4j) may also be applied. The extreme Reisner-Nordstrom geometry in which e = M is 
particularly interesting. In that case the two zeroes of /(r) at r± coincide and / vanishes quadrati- 
cally. The possible singularities of the stress-energy on the extreme Reissner-Nordstrom horizon are 
correspondingly more severe, behaving like s~^. Nevertheless a power series expansion near s = 
shows that it is possible to obtain a regular stress-energy on the horizon with ()3.36|) even in the 
extreme Reisner-Nordstrom case with appropriate choice of auxiliary field integration constants, in 
agreement with the direct numerical results of [4l|. Thus the anomaly induced effective action can 
reproduce the correct behaviors of the quantum stress tensor also in this case, in contrast to the 
approximation of 35|, 1361 . l37| which predicts a divergent stress-energy on the horizon of an extreme 
Reisner-Nordstrom black hole as the only possibility. As it is not possible to solve the auxiliary 
field equations analytically in the black hole geometries when e 7^ 0, and a numerical solution is 
required, we postpone a full discussion of these results to a future publication J49l |. 

Finally, this study the anomalous effective action, auxiliary fields, and stress-energy could be 
extended to rotating Kerr black hole geometries. The previous analytical approximations of 



se 
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43] do not apply to this case. However there is no problem of principle in computing the 
stress tensor ()3.36() in stationary spacetimes possessing only axial symmetry. Because of the lower 
symmetry, solutions of the auxiliary field equations (|3.1(ij) which are functions of both r and 6 
can be sought. Hence the solutions characterizing macroscopic stationary quantum states with 
axial symmetry are considerably richer than the Schwarzschild case. However, a simple spherically 
symmetric homogeneous solution of the second order equation D </:> = which generalizes the In / 
solution found in the the Schwarzschild and de Sitter cases is 

ip{r) ocln[ ^~^+ j , (5.19) 

for 



r± = M± Vm^ - a2 , (5.20) 

with a the the angular momentum per unit mass of the Kerr solution. Since (|5.19|) diverges 
logarithmically on the horizon at r = r^ , similar conformal behavior of the stress tensor as r ^ r+ 
is expected as that found in the Schwarzschild case, which clearly merits a separate investigation. 
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VI. CONCLUSIONS 

The suggestion that general relativity is a low energy effective field theory, comparable to the 
chiral EFT of pions in the strong interactions has been expressed previously 2] . In this paper we 
have extended the EFT framework from the weak field expansion around flat space considered in ^j 
to arbitrary curved spacetimes, provided only that the curvature invariants are small compared to 
the Planck scale. The EFT expansion in inverse powers of the high energy cutoff scale of order Mpi 
requires a systematic classification of relevant and irrelevant operators in the general coordinate 
invariant effective action. According to this classification, based on the behavior under global 
Weyl rescalings [lO| , the classical Einstein-Hilbert action should be supplemented with certain well- 
defined non-local terms associated with the trace anomaly of massless fields. These additional terms 
are non-local in terms of the metric, cf. (|3.12|) but nevertheless because they scale logarithmically 
with distance are relevant in the infrared. The anomalous terms are unsuppressed by any inverse 
power of the ultraviolet cutoff scale Mpi, and do not decouple for E <^ Mpi. Hence they are 
required in the Wilson effective action (|3.15|) for gravitational physics at macroscopic distances L ^ 
Lpi. This constitutes an non-trivial infrared modification of Einstein's theory which is completely 
consistent with the Equivalence Principle. 

The local form of the anomaly introduces two scalar auxiliary fields, (/? and ■0, into the EFT 
for low energy gravity, one for each of the non-trivial co-cycles of the Weyl group. These scalar 
fields satisfy the fourth order massless wave equations ()3.16() . and are only very weakly coupled to 
ordinary matter and radiation. In the flat space limit they decouple entirely. The auxiliary fields 
are new local scalar degrees of freedom of the gravitational field, not contained in classical general 
relativity. 

The anomalous terms in the effective action also generate two new conserved tensors (|3.38|) and 
(|3.39|) . local in the auxiliary fields, which contribute to the stress-energy tensor of low energy gravity 
H3.36() . as a source for Einstein's equations in vacuo. Since they arise from the long range effects of 
the fluctuations of massless fields, and require boundary conditions for their complete specification, 
these additional terms carry information about the global macroscopic effects of quantum matter. 
In the approximation where the auxiliary scalar fields are treated as classical potentials, they may 
be viewed as macroscopic order parameter fields of low energy gravity. 

Because of the conformal properties of spacetime horizons, the anomalous terms in the effective 
action and stress-energy become important in their vicinity, characterizing the allowed behaviors 
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of the stress tensor there. Near an horizon the most important terms in the effective action (|;-i.l5|) 
are those which scale positively in the ultraviolet short distance limit under conformal rescalings, 
while nevertheless keeping local curvature invariants fixed and small compared to the Planck scale. 
Only the Sanom term in (|3.15|) has this property, scaling logarithmically in both the ultraviolet 
and infrared limits. Because the anomalous terms depend non- locally on the geometry through 
the two auxiliary scalar fields 99 and ip, they can be significant near horizons even for very much 
sub-Plan ckian local curvatures, where the EFT description should continue to apply. Indeed al- 
though local curvature invariants remain small on horizons, the auxiliary scalar fields depend in 
general on coordinate invariant global quantities such as ln(— i^"i^a) which diverge on the horizon 
when the spacetime possesses a timelike Killing field K"^ which becomes null there. Whether or 
not this divergence actually occurs in the stress-energy of a given state is a matter of boundary 
conditions fixed globally over the entire spacetime. The existence of a topological current and as- 
sociated Noether charge are characteristic of these global effects of the quantum state. The generic 
singular behavior of the stress-energy on horizons is thereby characterized in terms of spacetime in- 
variant scalar order parameters in the low energy EFT, which is completely consistent with general 
coordinate invariance. 

In the (semi-) classical EFT framework where both the metric and auxiliary scalars are treated 
as classical fields, Eqs. (|3.16|) may be solved classically and provide a new method to calculate 
the approximate quantum stress tensor of conformal matter by classical means, bypassing the 
standard but cumbersome procedure of field quantization and renormalization of the quantum stress 
tensor. When states for which the stress-ener gy d oes not diverge on the horizon are considered. 



approximations such as those in 



35 



3e 



37l . l43l |. developed specifically for static backgrounds 
with a timelike Killing field yield a more accurate approximation to the expectation value of the 
stress-energy tensor. We conclude that in states with a finite stress-energy on the horizon the low 
energy Wilson effective action ()3.15() is not a particularly good approximation to the exact quantum 
effective action (|3.14|) . as indeed there is no reason why it should be in such states where quantum 
effects are negligibly small in any case. 

On the other hand when the scalar potentials diverge on the horizon, the associated terms in 
the stress-energy yield a good approximation to the vacuum expectation value of the stress-energy 
tensor. Compared to previous analytic approximations (|3.36|) with (|3.16j) has the considerable 
advantage of being completely general, and giving all possible allowed singular behaviors of the 
stress tensor on horizons. In this case the auxiliary scalar potentials may be viewed as horizon order 
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parameters. Even when the stress-energies diverge, they are characterized by auxihary fields with 
finite action, suggesting that such states should not be excluded a priori from physical consideration. 

We conclude that the effective action and stress-energy of the anomaly provides a consistent 
infrared modification of general relativity that incorporates long-range macroscopic quantum coher- 
ence and entanglement effects in a natural way. This low energy modification is actually required 
on general EFT grounds. The auxiliary fields of the anomalous action and stress-energy may be 
viewed as macroscopic condensates. 

The study of their effective stress-energy tensor in several special fixed spacetime backgrounds 
presented here clearly only begins a systematic treatment of their many possible other effects. In 
addition to extension to other backgrounds, particularly those with horizons, and a study of the 
scalar gravitational waves and long range effects these scalar potentials may generate, a rich variety 
of new solutions are to be expected when the Tab of Eqs. (|3.36|) is used as fully dynamical sources 
for Einstein's equations. Whether these effects can be used to test the theory against astrophysical 
observations is a most important issue for further investigation. In the case of divergent behavior of 
the stress-energy near the classical horizon a significant deviation of the geometry from the vacuum 
solutions of the purely classical Einstein equations is predicted. This raises the distinct possibility 
of qualitatively changing the global structure of both classical black hole and cosmological solutions 
|5fll |. The existence of an abelian conserved Noether charge in the EFT suggests comparison with 
a number current similar to that of non-relativistic many-body systems at a microscopic level {5l| . 
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APPENDIX A: STRESS-ENERGY IN STATIC, SPHERICAL, VACUUM SPACETIMES 

We give in this Appendix the explicit forms for the diagonal components of the two conserved 
tensors EJ' and F^*, cf. Eqns. (|3.38|1 and (|3.39j) . for static, spherically symmetric spacetimes with 
line elements of the form H4.17(l . The auxiliary fields are assumed to be of the form, 



If = if{r) + 2pKt , 
ip = ip{r) + 2p' Kt , 



(Ala) 
(Alb) 



with K the surface gravity at the horizon where / vanishes. All components of the two tensors 



were generated by MathTensor programs 
conservation and trace conditions (13. 4( 



521 ] ■ which were separately checked for satisfying the 



The time components of the two tensors so obtained are: 
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and 



Ft * = $(v^"V' + ^V" - v^'V") + ^ ( /' - - ) (v^'V + 9^V") + ^ 
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(A3) 



where primes denote differentaition with respect to r. 
The radial pressures of the two stress tensors are: 
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Finally the tangential pressures of the two tensors are: 
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The flux components E^^ and Fj'' are determined by the conservation Eqn., 
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(A8) 



to be proportional to 1/r^ with a proportionality constant given by the luminosity L as in (jS-lZf) 
of the text. For the general surface gravity k, T^*" = —L/A-Kr"^ with 



L = \Qttk. {bpq + bp'q + b'pq) . 
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<r/> 60 




FIG. la: The expectation value (Tj *) of a conformal scalar field in the Boulware state in Schwarzschild 
spacetime, as a function of s = ^—^f^ in units of 7r^r^/90. The solid curve is Eq. 13.36|1 with (|5.12b|) and 
c^ = — ^,d„ = ||, the dashed curve is the analytic approximation of [36|, and the points are the numerical 
results of jsj. 



2 


r\ 


1 


\ \ ^ 









/ : 


T/> -1 




1 

/ 
/ 


-2 




1 


-3 




1 

1 : 


-4 




' : 







10 



FIG. lb: The radial pressure (T^"") of a conformal scalar field in the Boulware state in Schwarzschild 
spacetime. The axes and solid and dashed curves and points are as in Fig. El 
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<V> 




FIG. Ic: The tangential pressure {Tg^) of a conformal scalar field in the Boulware state in Schwarzschild 
spacetinie. The axes and solid and dashed curves and points are as in Fig. ^| 
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FIG. 2a: The expectation value {T^ *) of a conformal scalar field in the Hartle-Hawking state in Schwarzschild 

spacetinie as a function of s = ^~~Kr' ™ units of 7r^T^/90. The solid curve is Eq. H3.36|l with Coo — 0.035, c^^ = 

-1.5144, p = 1.5144,^9 = q' = 2,d^ = 1.0262, d„ = 1.0096, p' = -1.0096^ the dashed curve is the analytic 

approximation of [43J and the data points are the numerical results of |40|. 
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FIG. 2b: The expectation value {Tf *) of a massless Dirac field in the Hartle-Hawking state in Schwarzschild 
spacetime. The sohd curve is Eq. (|3.36|1 with Cqo = 0.035, c^ = —1.571 = —p, q = q' ~2, d^o — 0.6059, d^ = 
0.6109 = — p', the dashed curve is the analytic approximation of |43l | and the data points are the numerical 
results of |41| . 
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FIG. 2c: The expectation value (Tj *) of the electromagnetic field in the Hartle-Hawking state in Schwarzschild 
spacetime. The sohd curve is eq. (|3.36|l with Coo = —0.0237, c„ = 1.0584 = p,q = q' = 2,doo = 
— 0.2716. d.. = 0.7644 = v'. and the dashed curve is the accurate analvtic fit of i42il to the numerical 
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FIG. 3a: The expectation value (Tj*) of a massless, conformal field in the Unruh state in Schwarzschild 

spacetinic. The solid curve is Eq. (I3.3()|l with Coo = 0,c^ = — 1.6500, p = 1.9192,(7 — q' — 2,doo — 

1.4441, d^ = 1.3244, p' = —1.0552, and the dashed curve is the polynomial approximation of 47] which is 

an accurate fit to the numerical results of [45|. 
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FIG. 3b: The expectation value (T^ *) of the electromagnetic field in the Unruh state in Schwarzschild 
spacetime. The solid curve is eq. ()3.36(l with Coo — 0,c^ — 0.5594, p — —0.4984,(7 — q' — 2,doo = 
— 0.2595, cf^ = 0.4125, p' = —.3514, and the dashed curve is the polynomial approximation of 48] which is 



